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FOREWORD 


This  is  the  final  engineering  report  of  the  Flight  Trajectory  Control  in¬ 
vestigation  program  under  Contract  F33615-78-C-3607  of  Air  Force 
Wright  Aeronautical  Laboratories,  Flight  Dynamics  Laboratory,  to  deter¬ 
mine  the  ,'easibility  of  applying  optimal  control  theory  to  four-dimension¬ 
al  guidance  of  a  transport  aircraft.  The  investigations,  in  conformance 
with  techniques  previously  developed  for  four-dimensional  guidance  of 
transport  aircraft  using  classical  control  system  theory,  included  appli¬ 
cation  of  optimal  algorithms  to  both  generation  of  a  three-dimensional 
reference  path  trajectory  through  a  sequence  of  waypoints  specified  in 
space  coordinates,  and  to  the  guidance  commands  that  controlled  the 
aircraft  along  the  reference  trajectory  at  a  velocity  necessary  to  achieve 
a  specified  time  of  arrival  at  the  waypoints . 

The  general  result  of  this  investigation  is  that  the  anticipated  degree  of 
improvement  in  aircraft  control,  using  optimal  rather  than  classical  algo¬ 
rithms,  could  not  justify  the  large  increase  in  computational  require¬ 
ments  of  the  on-board  avionics  computer  and  the  associated  increased 
costs  inherent  in  state-of-the-art  computer  hardware  technology. 

A  plan,  specifically  generated  as  part  of  this  investigation  to  demon¬ 
strate  the  optimal  algorithms  by  both  man-in-the-loop  simulation  using  a 
transport  fixed-base  cockpit  and  by  flight  test  in  a  transport  aircraft, 
can  be  implemented  for  demonstration  of  the  existing  classical  guidance 
algorithms . 

The  work  of  this  report  is  based  on  previous  engineering  development 
of  four-dimensional  guidance  of  transport  aircraft  using  classical  algo¬ 
rithms  as  reported  in  Technical  Report  AP’FDL-TR-77-120 . 

The  principal  contributors  to  this  report  are  the  project  engineer,  Dr. 
M.  Bird,  Dr.  R.  Fredricks,  and  Dr.  J.  Kamiarker,  Systems  Control, 
Inc.  (Vt).  Detailed  analyses  were  provided  by  Mr.  J.  Ring. 

The  investigations  of  this  report  were  accomplished  during  the  period 
May  1978  through  October  1980. 
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FINAL  REPORT 

FLIGHT  TRAJECTORY  CONTROL  INVESTIGATION 


INTRODUCTION 


This  report  details  the  results  of  the  Flight  Trajectory  Control  Invest¬ 
igation  (FTCI)  conducted  by  Lear  Siegler,  Instrument  Division,  Grand 
Rapids,  Michigan  under  Air  Force  contract  F33615-78-C  -3607. 

The  intent  of  the  FTCI  program  was  to  develop  improved  algorithms  in¬ 
corporating  optimal  control  techniques  for  guidance  and  control  of  a 
transport  aircraft.,  to  verify  these  algorithms  by  simulation,  and  to  de¬ 
sign  a  demonstration  flight  test. 

The  basic  algorithms,  developed  under  previous  Air  Force  contracts1, 
implemented  classical  control  techniques  for  real-time  generation  of  a 
flight  path  through  a  sequential  series  of  waypoints  defined  in  four 
dimensions  (x,  y,  z,  t)  and  for  control  of  a  transport  aircraft  to  the 
generated  flight  path. 

These  basic  algorithms  had  been  demonstrated,  using  fixed-base  cock¬ 
pit,  man-in-the-loop  simulations,  t.o  Drovide  excellent  functional  and 
operational  guidance  and  control  capability  for  a  transport  aircraft. 
However,  these  algorithms  were  not  'optimal',  in  the  sense  of  minimizing 
a  specified  cost  function  that  could  be  defined  in  terms  of  spacte  posi¬ 
tion  and  time  control  performance  trade-offs. 

A  principal  issue  in  this  FTCI  program  was  the  cost  effectiveness  of 
any  optimal  control  algorithms ,  as  compared  to  the  classical  control 
algorithms.  The  increased  cost  of  applying  optimal  control,  as  mea¬ 
sured  by  increased  processing  speed  and  memory  requirements  of  the 
airborne  avionics  computer  to  handle  the  more  sophisticated  algorithms, 
must  be  paid  by  increased  performance  as  measured  by  deviations  of 
the  actual  aircraft  flight  path  from  the  specified  waypoints  that  define 
the  desired  flight  path , 

Early  in  the  program  it  became  apparent  that  the  cost  for  initially  de¬ 
veloping  the  optimal  algorithms  for  generation  of  the  trajectory  in  four 
dimensions  would  greatly  exceed  the  contractual  resources.  The  pro¬ 
gram  was  redirected  to  develop  the  optimal  algorithms  to  control  the  air¬ 
craft  to  a  four-dimensional  trajectory  defined  by  the  existing  classical 
algorithms.  The  initial  control  algorithms  developed  in  the  form  of  a 


feasibility  Study  for  Integrated  Flight  Trajectory  Control  (Airlift). 
LSI  Report  #  GRR- 008-01 77A ,  AFFDL-TR-77-120 ,  contract  F33615-74-C- 
3083. 
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linear  quadratic  estimator,  did  not  show  any  improvement  in  control  per¬ 
formance  when  compared  to  the  classical  algorithms.  As  a  result,  tne 
flight  test  plan,  while  useful  for  demonstrating  either  optimal  or  classi¬ 
cal  algorithms,  is  based  on  an  avionics  system  mechanization  that  incor¬ 
porates  only  the  classical  algorithms. 

Although  the  major  findings  of  this  program  indicate  that  optimal  control 
applied  to  the  basic  Integrated  Flight  Trajectory  Control  functions  is 
not  yet  within  the  state-of-the-art,  continuing  advances  in  methods  of 
implementing  optimal  algorithms  for  aircraft  guidance  and  control  and  in 
airborne  processor  mechanizations  may  yet  moke  the  FTCI  concepts  not 
only  feasible  but  desirable. 


PROGRAM  SUMMARY 


2.1  PROGRAM  PLAN 

Two  ksy  elements  of  the  general  plan  to  achieve  the  tasks  of  the 
FTCI  program  were  (1)  to  utilize  the  assets  accruing  from  the  preced¬ 
ing  IFTC  program  program  and  (2)  to  obtain  specialized  technical  sup¬ 
port,  through  a  subcontract,  from  System  Control  Incorporated  (SCI)  of 
Florida. 

The  major  program  tasks  included  development  of  optimal  trajectory  gen¬ 
erator  and  aircraft  control  algorithms,  verification  of  these  algorithms 
by  man-in-the-loop  simulation  using  a  fixed-base  cockpit,  and  design  of 
a  flight  test  program  to  demonstrate  the  performance  of  the  algorithms 
in  a  typical  operational  system 

The  assets  or  background  for  this  program  that  were  derived  from 
IFTC  included  classical  algorithms  that  were  the  functional,  but  not.  ne¬ 
cessarily  the  performance,  equivalent  of  the  desired  optimal  algorithms. 
In  addition,  the  assets  included  a  fixed-base  transport  cockpit  simulator- 
equipped  with  man-machine  interface  specifically  developed  for  the  oper¬ 
ational  requirements  of  IFTC  mission  management  systems. 

The  subcontract  with  SCI  (Vt)  was  designed  to  provide  Lear  Siegler, 
Inc.  (LSI)  with  experienced  technical  advice  and  support  in  the  appli¬ 
cation  of  nonlinear  optimal  control  theory  to  the  general  task  of  algor¬ 
ithm  development  and  to  provide  novel  concepts ,  developed  by  SCI ,  for 
threat  modeling  and  avoidance  capabilities  for  the  trajectory  generator . 
The  expertise  of  SCI  (Vt)  in  these  areas  is  based  on  successful  related 
programs.  A  typical  example  is  the  Advanced  Weapons  Management 
Systems  (AWMS)  for  the  Navy.  In  this  program  SCI  (Vt)  developed  a 
dynamic  programming  solution  for  determining  at  optimum  flight  path 
and  allocation  of  other  resources  for  tnreat  avoidance  by  a  manned  air¬ 
craft  on  any  enemy  defense  penetration  mission. 

2.2  PROGRAM  HISTORY 

The  initial  task  in  developing  the  algorithms,  as  shown  in 
Figure  1,  was  to  develop  a  first-cut  approach  to  applying  optimal  con¬ 
trol  theory  to  the  trajectory  generator. 

The  approach,  described  in  3.1,  was  straightforward  end  amounted 
basically  to  adapting  the  classical  two-point  boundary  value  formulation 
to  a  six-degree-of-freedom  point  mass  aircraft  with  all  the  structural, 
maneuver,  and  mission  constraints  accounted  for  by  using  penalty  func¬ 
tions  (both  terminal  and  enroute)  in  J,  the  overall  cost  function.  An 
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eight- dimensional  discrete  state  vector  (x)  was  formulated  consisting  of 
east,  north,  and  vertical  aircraft  displacements,  true  airspeed,  flight 
path  angle,  heading  angle,  vehicle  mass,  and  elapsed  time.  The  last 
state  variable  was  introduced  because  ground  track  distance  was  used 

as  the  independent  variable  rather  than  time.  The  control  vector  (u) 
chosen  consisted  of  three  elements,  namely,  angle  of  attack  (alterna¬ 
tively  pitch),  craft  roll  angle,  and  effective  engine  power  setting  or 
EPR. 

A  zero-th  order  solution  for  the  complete  trajectory  including  a  terminal 
value  (R  for  the  ground-track  variable  (r)  was  assumed  to  be  obtain¬ 
able  from  the  existing  conventional  IFTC  trajectory  generator  which  has 
already  been  developed.  The  solution  for  the  state  (x),  costate  (X) 

(also  eight- dimentional) ,  and  control  vectors  over  the  complete  trajec¬ 
tory  was  obtained  using  a  steepest  descent  technique  to  minimize  the 
gradient  of  the  Hamiltonian  (H)  at  each  discrete  state  of  each  iteration. 
The  control  vector  trajectory  was  continually  perturbed  via  the  steepest 
descent  algorithm  until  V  H,  the  gradient  of  H  with  respect  to  u,  be¬ 
came  sufficiently  small  at  each  of  the  assumed  N  discrete  track  stages 
in  the  profile  that  a  stopping  criteria  was  satisfied.  This  process 

yielded  a  candidate  optimal  solution  [x{r} ,A{r} ,u{r] ]  for  the  trajectory, 
valid  from  r=0  to  r=RQ.  New  candidate  optimal  solutions  were  obtained 
using  an  outer  iteration  loop  where  Rq  was  perturbed  from  the  conven¬ 
tional  IFTC  generator  value.  By  calculating  the  actual  cost  associated 
with  each  candidate  optimal  solution,  the  true  optimum  could  be  deter¬ 
mined. 

It  was  felt  that  this  algorithm  was  amenable  to  a  real-time  airborne 
implementation  due  to  the  speed  and  small  size  of  current  military  micro¬ 
processors.  However,  an  algorithm  review,  by  subcontractor  SCI,  indi¬ 
cated  that  the  task  of  implementing  this  algorithm  in  our  simulator 
would  exceed  substantially  the  manpower  and  dollar  resources  aHotted  to 
the  FTCI  program.  In  addition  to  excessive  algorithm  development 
costs,  SCI  indicated  that  the  penalty  function  approach  to  constraint 
satisfaction  was  too  "soft",  even  though  it  did  make  the  full  two-point 
boundary  value  quite  suitable  for  the  steepest  descent  solution.  Al¬ 
though  the  constraints  could  be  forced  to  be  satisfied  more  exactly  if 
more  complex  optimization  methods ,  such  as  gradient  projection ,  were 
used  to  augment  the  steepest  descent  approach,  proper  investigations  of 
these  techniques  would  only  exacerbate  the  time,  manpower,  and  financ¬ 
ial  overload  of  the  program. 


As  an  alternative,  SCI  suggested  implementing  a  piecemeal  procedure 
that  would  generate  s_,boptimal  flight  profiles.  This  approach  would 
first  determine  a  horizontal  plane  profile  which  avoided  all  threats  and 
satisfied  all  the  waypoint  and  endpoint  constraints.  This  profile  could 
then  be  displayed  to  the  pilot  for  his  approval.  Next  would  be  tu  use 
an  energy-state  optimal  control  algorithm,  similar  to  that  being  pursued 
by  Erzberger .  et  al.2,  to  calculate  the  vertical  profile,  speed,  and  time 
schedules  associated  with  the  accepted  horizontal  trajectory.  This  ap¬ 
proach  would  take  maximum  advantage  of  the  threat  avoidance  (minimum 
exposure)  horizontal  profile  generator  previously  developed  during  the 
AWMS  program  and  of  the  considerable  previous  work  based  upon 
Erzberger's  energy-state  algorithm,  as  reported  in  the  literature. 

The  piecemeal  approach,  however,  had  the  disadvantage  of  limited  in¬ 
herent  performance  improvement  capability. 

The  problem  was  resolved  by  rescop  mg  the  FTCI  program  to  concen¬ 
trate  the  remaining  resources  in  developing  the  linear  quadratic  optimal 
control  algorithms  for  control  of  the  aircraft  to  a  trajectory  developed 
using  the  IFTC  classical  techniques  but  modified  with  increased  surviv¬ 
ability  capability  by  incorporating  algorithms  developed  by  SCI.  These 
algorithms  are  defined  in  Section  3.3. 

The  development  of  the  linear  quadratic  optimal  controller,  as  detailed 
in  Section  3.2,  began  with  definition  of  dimentionless  force  and  moment 
equations  for  the  KC-135  aerodynamics.  These  equations  were  derived 
by  simplification  of  the  standard  aerodynamic  equations  defined  by  Boe¬ 
ing  Military  Aircraft  Division  for  the  KC-135  aircraft. 

With  these  aerodynamic  forces  and  moments  the  equations  of  motion  in 
six  degrees  of  freedom  were  generated  for  simulation  of  the  KC-135  air¬ 
craft.  The  state  model  of  the  aircraft,  for  this  purpose,  was  defined 
in  terms  of  nine  tate  variables:  airspeed,  angle  of  attack,  angle  of 
sideslip,  and  boay  axis  angular  rates  in  inertial  space  and  attitude. 
This  aircraft  model  was  simulated  by  hybrid  comp’  ter  and  the  simulation 
was  verified  by  comparing  the  response  of  the  simulated  system  to  simu¬ 
lated  aileron  and  elevator  deflections  to  calculated  responses  of  identical 
analytic  inputs. 

The  inner  control  loop  of  the  aircraft  simulation  was  closed  by  design 
and  simulation  of  the  stability  augmentation  system  (SAS)  using  attitude 
and  rate  feedback. 


2H.  Erzberger,  et  al.,  "Fixed  Range  Optimum  Trajectories  for 
Short-Haul  Aircraft",  NASA  TN-D-8115,  December  1975. 
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The  SAS  and  autothrottle  design  were  based  on  the  current  digital 
flight  control  system  being  installed  on  the  Speckled  Trout  aircraft.. 

To  complement  these  simulations  an  IR&D  project  rehosted  the  1FTC  al¬ 
gorithms  from  our  IBM  370  to  a  PDP-11/70  laboratory  computer.  This 
improved  operating  time  because  the  IBM-370  was  entered  through  an 
often  unreliable  time-sharing  real-time  monitor  system  while  the  PDP- 
11/70  is  a  dedicated  machine  for  avionics  systems  simulation.  In  addi¬ 
tion,  the  controller/display  unit  of  the  cockpit  simulator  was  modified  to 
provide  improved  operability. 

To  reduce  the  complexity  of  the  linear  quadratic  (LQ)  controller,  the 
aircraft  equations  of  motion  were  reduced  to  point  mass  rather  than  the 
six-degree-of-freedom  equations  employed  for  the  simulation.  These 
point  mass  equations  were  linearized  and  used  to  develop  a  simplified 
form  of  the  aircraft  state  equations.  To  analytically  determine  the 
steady-state  performance  of  the  optimal  controller  to  a  deterministic  dis¬ 
turbance  such  as  a  change  in  horizontal  wind  magnitude  and  direciton, 
the  simplified  aircraft  state  model  was  developed  as  a  function  of  the 
following  state  variables,  airspeed,  heading,  flight  path  angle,  position 
in  along-/across-track  coordinates,  and  altitude.  The  driving  functions 
were  modeled  as  assumed  rates  of  change  of  wind  magnitude  and  direc¬ 
tion.  For  this  input,  an  analytic  solution  of  the  state  perturbation  ma¬ 
trix  and  control  correction  matrix  were  obtained  to  demonstrate  feasi¬ 
bility  of  the  LQ  optimal  controller,  with  a  constant  value  for  the  state 
(Q)  and  control  (R)  gain  feedback  weighting  matrices. 

a.  A  diagonal  form  of  the  weighting  matrices  was  selected  to 
allow  independent  selection  of  values  for  individual  states  and 
controls . 

b.  Initial  values,  Qq  and  Rq,  were  established  for  a  particular 

flight  condition  by  assigning  maximum  values  to  the  plant 
state  and  computing  maximum  allowable  controls  corresponding 
to  these  states. 

c.  The  components  of  the  matrices  were  varied  to  determine  the 
effect  on  the  closed-loop  system  optimal  controller  stability  by 
eigenvalue  analysis . 

d.  A  sensitivity  analysis  determined  that  three  elements  of  the  Q 

matrix  (6V_  ,  ,  6tb _ and  6 X  )  which  penalized  errors  in 

max  max  max  r 

velocity,  heading,  and  flight  path  angle  must  be  modified  as  a 

function  of  flight  profile  to  achieve  the  desired  response. 


e.  The  weighting  matrices  were  scheduled  as  a  lunation  of  flight 
condition  by  selecting  a  performance  index  for  a  sample  of 
representative  flight  conditions  that  resulted  in  a  critically 
damped  closed-loop  system. 

The  simplification  of  the  controller  by  incorporating  a  point  mass  air¬ 
craft  model  and  a  stability  augmentation  control,  rather  than  an  attitude 
flight  control  system,  minimised  computation  time  in  the  airborne  proces¬ 
sor  but  also  degraded  probable  performance.  The  cross-coupling  terms 
derived  by  the  controller  would  make  aircraft  transient  responses  to  pi¬ 
lot  overrides  and  to  scheduled  gain  changes  as  a  funciton  of  flight  con¬ 
dition  unacceptable  to  the  pilot.  No  further  funds  were  expended  dir¬ 
ectly  for  optimal  control. 

A  flight  test  program  to  mechanize  the  Speckled  Trout  test  aircraft  for 
demonstration  of  IFTC  characteristics  of  time  control  and  in-flight  mis¬ 
sion  redirect  was  generated.  This  flight  test  plan  was  designed  as  a 
three-phase  program  because  of  the  limited  modification  time  availability 
of  the  Speckled  Trout  test  aircraft  and  because  of  the  difficulty  in  com¬ 
plete  single  step  replacement  of  large  sections  of  the  avionics  on  an 
operational  aircraft.  This  flight  demonstration  plan,  as  described  in 
Section  4,  includes  the  avionics  equipment  description,  the  flight  test 
scenarios,  and  the  installation  plan. 
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3. 


ALGORITHM  DEVELOPMENT 


Algorithm  development  was  initiated  by  a  preliminary  definition  of  a  pre¬ 
ferred  approach  to  the  solution  based  on  the  concepts  of  the  original 
proposal.  This  initial  approach  was  expanded  in  further  detailed  de¬ 
velopmental  analyses  for  both  the  optimal  trajectory  generator  and  opti¬ 
mal  controller. 

The  optimal  control  problem  was  based  on  the  classical  algorithms  deve¬ 
loped  on  the  previous  Flight  Trajectory  Control  System  Investigation. 
These  algorithms  defined  a  flight  path  of  an  aircraft  which  started  from 
an  "initial  point"  and  terminated  at  a  desired  "end  point".  A  "point"  is 
defined  as  a  four-dimensional  (space,  time)  coordinate  set  that  specifies 
a  trajectory  as  well  as  the  state  vector  of  the  aircraft  (heading,  velo¬ 
city,  and  flight  path  angle)  with  respect  to  the  specified  trajectory. 

The  objective  of  this  present  Flight  Trajectory  Control  Investigation  was 
to  develop  algorithms,  based  on  optimal  control  techniques,  that  would 
most  nearly  control  the  aircraft  to  the  optimized  trajectory.  Because 
the  constraints  on  optimization  were  the  aircraft  dynamics,  the  optimized 
trajectory  had  to  most  nearly  conform  to  the  capabilities  of  the  aircraft 
while  achieving  the  desired  mission  requirements;  and  the  optimized  con¬ 
troller  had  to  minimize  the  actual  excursions  of  the  aircraft  from  the  tra¬ 
jectory  during  flight. 

Classical  deterministic  optimal  control  algorithms  were  developed  by  ap¬ 
plying  the  principles  of  variational  calculus  to  the  initial,  non-optimal, 
trajectory  generator  and  control  algorithms  developed  on  the  IFTC  pro¬ 
gram.  The  high  computational  workload  required  to  implement  these 
classical  optimal  control  techniques  has,  to  date,  made  them  unsuitable 
in  applicability  to  airborne  processors,  The  rate  of  advance  in  speed 
of  digital  processors,  however,  makes  these  optimal  algorithms  feasible 
for  advanced  avionics  system  mechanizations. 

3.1  TRAJECTORY  GENERATOR  DEVELOPMENT 

3.1.1  Development  Method 

A  six-degree-of-freedom  point  mass  model  is  generated  for  the 
aircraft  equations  of  motion.  Thus  x,  y,  z,  v,  y  and  tj)  are  automatic¬ 
ally  state  variables  with  x  east,  y  north,  z  up,  v  the  true  airspeed,  y 
the  flight  path  angle  and  4/  the  heading.  To  these  add  aircraft  mass, 
m,  and  since  ground-track  distance,  r,  is  the  independent  variables, 
time,  t,  is  also  a  state  variable. 
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The  control  vector,  u,  consists  of  three  elements:  angle  of  attack,  a, 
roll  angle,  <j>,  and  effective  power  setting,  n.  The  latter  control  vari¬ 
able  is  a  number  between  zero  and  100  per  cent  and  represents  the  ro¬ 
tor  RPM  (effective  rotor  RPM  for  multi-engine  craft)  normalized  to  the 
maximum  allowable  RPM.  Given  a  power  setting  command,  n,  the  auto¬ 
pilot  can  readily  set  the  throttle  levers  to  achieve  the  desired  RPM. 

The  cost  function  (J)  includes  terminal  costs  (K)  and  integral  (path) 
losses  (I).  The  integral  losses  incorporate  constraints  such  as  maxi¬ 
mum  and  minimum  velocities  or  altitudes,  and  avoidance  regions  such  as 
cylindrical  volumes  about  some  given  ground  coordinate.  Penalty  func¬ 
tions  incorporate  all  such  path  constraints  into  J. 

The  state  (x)  and  costate  (X)  vector  difference  equations  are  solved 
iteratively  for  each  path  increment  Ar,  where  the  magnitude  of  Ar  will 
be  determined  by  simulation.  A  probable  "safe"  value  of  Ar  suitable  for 
transport  dynamics  is  '>-0.1  nautical  mile.  The  number  of  the  iterative 

o 

multipass  stages  in  the  solution  of  the  algorithm  is  designated  by  N^. 

This  number  is  fixed  for  a  given  £  independent  of  k.  The  nominal 
initial  trajectory  has  a  total  ground  track,  RQ.  Thus  the  number  of  all 

iterations  of  the  first  pass  of  the  algorithm  is 


where  Ar  is  small  enough  to  give  accurate  results  and  selected  to  make 
Nl  integer. 

The  index  (i)  on  discrete  increments  of  ground  track  then  runs  from 
i=0  to  i=N  for  the  kth  iteration  of  the  £th  pass.  The  state,  costate, 
and  Hamiltonian,  are  solved  for  at  each  r..  The  control  vector,  how¬ 
ever,  is  considered  to  vary  more  slowly  and  the  steps  in  the  discrete 
control  vector  occur  only  ever  j=10i  units  in  discrete  ground  track  with 

the  dummy  variable  i  running  from  0  through  — 

The  relatively  simple  steepest  descent  algorithm  for  updating  the  control 

vector  (u.  .  )  at  the  (k+l)st  iteration  of  a  particular  pass  was  imple- 
J  >  k 

mented.  Conjugate  gradient  and  gradient  projection  techniques,  al¬ 
though  improving  convergence  time,  are  too  complex  for  incorporation  in 
an  airborne  computer. 
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In  addition  to  updating  the  control  vector  (u.  ,  )  by  a  steepest  descent 

J  > 

iterative  approach,  an  ad-hoc,  multi-pass  technique  will  be  used  to 

o 

slightly  modify  N  ,  the  number  of  Ar  stages.  The  approach  of  replac¬ 
ing  time  as  the  independent  variable  with  ground  track  (r)  results  in  a 
problem  with  length  of  ground  track.  Thus  for  the  £th  pass,  a  change 
£ 

in  N  represents  an  alteration  of  Rq.  The  constraint  on  time  is  repre¬ 
sented  in  the  terminal  cost  function  (K)  as  a  term  cf  the  form 
^N^’k'V2  where  t^  is  the  desired  terminal  time,  t^£,k  is  the  value 
of  the  state  variable  "t"  at  the  terminal  state  N£  for  the  kth  iteration  of 
the  £th  pass,  and  ktt  is  some  non-negative  penalty  function. 

Then  the  complete  algorithm  requires  iteration  on  u.  .  0  at  each  stage  j 

for  a  given  pass  (£=Constant)  with  k  representing  the  iteration  number. 
When  a  stopping  criteria  is  satisfied,  the  iterations  for  terminated  at 

k=k*  and  the  candidate  optimal  control  discrete  sequence  { ^ }  is 

used  to  evaluate  an  overall  cost  function  J.  Next,  £  is  incremented, 

i.e,  ,  N  is  perturbed  (both  above  and  below  N1)  and  new  control  se¬ 
quences  calculated  via  iterative  steepest  descent  techniques. 

The  final  optimal  trajectory  is  then  taken  as  the  state  vector  generated 

by  the  optimal  choice  of  both  N  and  {u.  Letting  £*  denote  the 

o#  J » K  ^ 

optimal  terminal  distance  (via  R*=N  Ar),  the  control  sequence  for  the 
various  stages  is  given  by  {Uj  £*}  where  j=10i  and  i  =  0,1,2 . 

a* 

(N-JlLO) 

^  10  1 

£ 

The  optimal  solution,  N  ,  is  the  point  of  minimum  J,  the  cost.  Thus  J 

£ 

is  a  function  of  both  the  iteration  k  and  the  number  of  stages  N  ,  i.e., 
J=J  { k ,  Jd } . 

Only  a  limited  search  was  made  in  £,  i.e.,  £=1,2,3  was  used  when  N1 
was  as  given  previously,  N2=N1(1.1)  and  N3=N1(0.9)  say.  The  best  of 
these  three  was  used  to  define  the  final  optimal  flight  path. 


3.1.2  Dynamics  of  Problem 

The  aircraft  equations  of  motion  assuming  zero  sideslip  angle  and 
neglecting  round,  rotating  earth  effects,  can  be  written  as 


X 

Vcos{y  }sin{t|/}  +  w»in{a} 

(1) 

y 

= 

Vcos{y}cos  {iJj}  +  wcos{a} 

(2) 

• 

z 

= 

Vsin{y} 

(3) 

• 

V 

- 

^  cos  {a  +  CD  m  “  8sin{vl 

(4) 

V 

= 

[j  si“!“  +  V  +  2V  (V>‘  CL  i]c0,{lK 

(5) 

-  ^  cos{y} 

4> 

= 

[m  Sl“>“  *  “T1  *  2  tV)2  CL  i] 

(6) 

• 

m 

s 

U){T,V,Z} 

(7) 

The  ground  track  continuous  dynamics  may  be  derived  from  these  equa¬ 
tions  by  substitution  of  equation  8  for  the  time  differential,  dt.  Then 
equations  1  through  7  become  equations  9  through  15. 


Tfr  1 1!0)  - 0 


•  t  i  i  /W\  sin  of 

“"t*1  T  cv)  iSlfoj 


,  .  ,  .  wcos (O : 


tan{y} 


(12) 


dv 

dr 


dl 

dr 


d»Jj 

dr 


dm 

dr 


In  these  equations 


V 

Y 

4> 


m 


T{V,z,n}cos{a+aT} 
mVcos {y} 

_  gtan{y} 

V 

T{a.V,n}sin{o(+aT} 
mV^cos fyj 

*  cos  { <t»}  - 

T{z,V,rt}sin{a+aT} 

mVcos2{y} 

•  sin{4>} 

-m{z,V,n} 

Vcos {y} 


p{z}CD{z,V,a}S 

2rcos{y} 


p{z}CL{z,V,a}S 
2mcos { y } 


p{z}VCL{z,V,a}S 

2mcos^I'Y] 


(13) 


(14) 


(13) 


airspeed 

angle  of  velocity  vector  with  respect  to 
the  x-y  plane;  positive  y  corresponds 
to  a  climb  angle 

velocity  vector  heading  measured  in  the 
x-y  plane  clockwise  from  North  to  the 
projection  of  the  velocity  vector  in  the 
x-y  plane 

mass  of  aircraft 


a  =  angle  of  attack;  positive  a  corresponds 

to  aircraft  longitudinal  axis  pitched 
above  the  velocity  vector 

4>  =  roll  angle,  positive  right  wing  down 

T  =  aircraft  thrust 


a™  =  constant  angle  of  thrust  with  respect  to 
the  lift  line  of  the  aircraft 
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w,o 

wind  magnitude  and  heading  with 
respect  to  North 

P 

= 

atmospheric  density 

CD 

= 

drag  coefficient 

S 

= 

vehicle  reference  area 

e 

= 

gravitational  acceleration  (positive) 

The  complex 
atmospheric 

functional 

parameters 

interrelationships  of  the  KC-135  vehicle  and 
are  summarized  by  the  following  equations : 

CL 

= 

CL{a,  z,  //M,  q} 

(16) 

- 

ko  +  ki  a 

CD 

- 

CD{Cl’  m] 

(17) 

= 

CDq{#M}  +  k2 {y/M}  (CL  -  0.1) 2 

(18) 

T 

A 

NRT  •  TS 

(19) 

where 

#M 

= 

Mach  number 

q 

= 

dynamic  pressure 

NET 

— 

normal  rated  thrust  from  all  aircraft 
engines 

TS 

= 

effective  throttle  setting 

T 

= 

thrust 

In  equations 

12  through 

15, 

however,  thrust  (T)  and  fuel  flow  (uj)  are 

expressed  directly  in  terms  of  n,  the  effective  power  setting  that  de¬ 
notes  actual  rotor  RPM,  and  pressure  and  temperature  factors,  6  and  0, 
respectively .  Then 
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H 


'U  r> 


with 


T  =  #M} 

uj  =  6 Ve  f F{nVe,  #m} 


(20a) 


(20b) 


6  =  (X{z}/T0)  (1  +  0.2  m2)  (20c) 

6  =  (P{z}/P0)(l  +  0.2  (/M2)3*5  (20d) 


where  we  have  assumed  the  specific  heat  ratio  of  air  at  constant  pres¬ 
sure  to  that  of  air  at  constant  volume  to  be  1.4  and  where  t{z)/t0  and 
P{z}/P0  are  atmospheric  temperature  and  pressure  ratios  obtained  from 
the  1962  standard  atmospheric  model.  Measurement  of  the  stagnation 
pressure  and  static  pressure  allows  direct  determination  of  Mach  num¬ 
ber.  The  latter  is,  of  course,  linearly  relatable  to  V,  if  the  free 
stream  static  air  temperature  is  known. 


The  control  vector  (u)  is  three-element  vector  given  by 


U1  {t} 

a{t] 

u2  ft] 

- 

4»{t} 

u3{t} 

T{t} 

(21) 


In  the  preceding  equations  we  have  assumed  a  flat  earth  coordinate  sys¬ 
tem  with  x  east,  y  north,  and  z  vertical. 


Given  N  and  Ar,  the  discrete  versions  of  equations  8  through  15  are: 


fci+l 


t. 

i 


+ 


Ar 

vTcosTy~T 


(i=0,l , . . . , [ N— 1 ]  and  to=0) 


(8 1 ) 


x 


i+1 


x .  +  sin  {tin } 


w. sin{o. } 

+  V.  cosfy!} 

l  1  ’ 


(9') 


(i=0,l,.. . ,[N-1J) 


3.1.3  Cost  Function 


Assuming  an  N  stage  process  the  total  cost  function  (J)  is  the 
sum  of  the  terminal  (K)  and  integral  (I)  (summation)  costs.  We  write 

J  =  K  +  i 

where  K  is  defined  by 

K  =  (xN  -  x£p)  [KjjKxjj  -  x£p) 

and 

“  lep 

XN  "  XEP 
yN  "  ZEP 
ZN  '  ZEP 
VN  "  VEP 
YN  ‘  YEP 
%  ’  ^EP 

“n  "  "'ep 

with  the  subscript  "EP"  denoting  the  end  point. 

[K^j]  itself  is  an  8x8  diagonal  matrix  representing  the  terminal  cost 
weights . 


(XK  -  *Ep) 


[K] 


tt 


XX 


yy 


k  o 

zz 


w 


YY 


«w» 


K 


mm 


where  some  elements  such  as  K  or  K  may  be  zero,  reflecting  "don't 
care"  conditions .  1  YY 
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(22) 


(23) 


(24) 


(25) 


+ 


(26) 


I 


where 

Lz. 

1 

LV. 

1 

La. 

l 

A  possible  formulation  of 


N-l 

~a  ~ 

1 

"a 

=  l 

< 

[R] 

♦I 

o 

II 

71  . 

i_ 

n. 

i- 

+  Lz .  +  LV .  +  La . 

ill 


=  altitude  constraints 

=  velocity  constraints 

=  avoidance  constraints 

Lz.  is 


,ui{z  .  ,V . ,n . } 
fuel  i  i  i 

V.cos{y. } 
i  *i 


Lz . 


i 


k  exp{- (z . -zmin. )/o  }  +  k  ' 

•y  r  '  '  i  '  rr  *  7 


exp{+(z . -zmax . )/o  } 
r  1  1  z 

(27) 


where  a  =  200  feet  and  zmin.,  zmax.  are  determined  in  an  ad  hoc  man- 
z  11 

ner  from  the  initial  trajectory  solution.  Without  any  ATC  constraints, 
zmin.  might  be  set  equal  to  1000  feet  above  the  highest  obstacle  any- 

wnere  within  a  0.5°  forward  look  angle  of  the  aircraft.  Similarly,  with¬ 
out  any  ATC  constraints  zmax.  might  be  set  at  40,000  feet  MSL  for  all 

i.  Under  Air  Traffic  Control,  zmin.  might  be  set  at  a  flat  18,000  feet 

MSL  and  zmax.  v  equal  40,000  feet  MSL.  During  takeoff  and  approach 

segments,  zmin.  must,  of  course,  be  suitably  reduced.  The  demarcation 

boundaries  separting  such  segments  from  those  of  cruise,  cruise/climb, 
or  cruise/  descent  can  be  determined  in  terms  of  r.,  the  independent 

variable,  from  the  nominal  trajectory,  {x.,q;u.,q}. 


A  possible  formulation  of  LV.  is 


LV. 


ky  exp  - 


r (v.-v  .  )i 

i  min 

+  ky'  exp  + 

"(V.-V  )“| 

i  max 

a,, 

O,, 

L  v  J 

L  V  J 

(28) 
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. . -  £  ika: ; 


where 


V  .  =  1.2Vk  „ 

min  stall 

V  =  max  dash  velocity 

max 

av  =  5  kts 

A  possible  formulation  for  La.  is 

k 

_ xaya _ 

Ui  =  |-(x.-xa)“  +  , 

expL — v  J 

where,  typically,  a  might  be  50,000  feet  to  500,000  feet  and  the 

xy 

avoidance  cylinder  is  centered  at  (xa>ya)‘  These  might  be  the  x,y  co¬ 
ordinate  of  the  SAM  site  or  of  n  mountain  peak.  Similar  cost  functions 
may  be  formulated  for  other  avoidance  geometries. 

The  cost  function  for  fuel  (kfuel>  in  I  represents  the  penalty  for  using 

too  much  fuel  per  unit  of  distance.  The  symmetric  positive  semi- 
definite  matrix  R  represents  the  penalty  on  too  much  control  action.  It 
might  be  more  meaningful  to  penalize  control  deviations  from  the  original 
nominal  control  stragey,  i.e.  ,  use 


a .  -  a .  _ 

"  ail 

i  i,0 

1 

♦i  '  *  i , 0 

in  place  of 

•H 

-o- 

n .  -  n .  n 

n. 

i  i,0 

L 

where  it  is  assumed  the  nominal  control  (a,  q ; ^i , 0 ; ni , 0^  is  slowly  vary 

ing  vector  with  respect  to  path  length.  This  control  strategy  penalizes 
execessive  perturbation  that  requires  continuous  new  inputs  to  the  auto¬ 
pilot/autothrottle  . 

3.1,4  The  Discrete  Stage  Optimal  Control  Algorithm 

The  general  form  of  the  discrete  state  dynamic  equations  8' 
through  15'  may  be  written, 


(29) 
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(30) 


x0  given 


(i“0 , 1 » • • • [N-l ] ) 

For  a  given  pass,  N  is  considered  fixed.  The  Hamiltonian  (H.)  is  de¬ 
fined  as  a  scalar  quantity  given  by 


H. 

1 


I .  +  K.  . ,  A  {x . ,u . } 
l  l+l  l’  l 


(31) 


where  1.  is  an  element  of  the  integral  cost  function  as  defined  in  equa¬ 
tions  26  through  29,  and 


is  the  final  value  of  the  costate  vector 
satisfies  the  equation: 


[kn]  [xn“xep] 

In  general  the  costate  vector 


(32) 


9H 


i 


where  x  is  an  eight-state  vector. 

The  necessary  condition  for  a  maximum  of  H  is  that 


9H 


au 


a 


J  1=J 


where  u  is  a  three-element  control  vector,  given  in  equation  21,  which 
drives  the  state  vector. 

In  equation  34,  index  j  is  a  multiple  of  a  base  number  of,  say,  10, 
i.e.  , 


(33) 


(34) 


(j  =  0,  10,  20,  ...,  (N-10) ) 

This  relationship  of  j  to  i  is  selected  so  that  the  control  will  vary  only 
1/10  as  rapidly  as  the  state  or  co-state  for  a  given  Ar.  This  assump¬ 
tion  reduces  the  number  of  individual  necessary  equations  that  must  be 

solved  from  3N  to  3(y^). 
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4  f*J 

Our  total  number  of  unknowns  (given  x^  and  X^)  is  nN+nN+m  yy  where 
n=9  states  and  m=3  controls.  The  state  dynamics  given  by  vector  equa¬ 
tion  30  yields  nN  scalar  equations  as  does  vector  equation  33.  Finally, 

N 

the  necessary  vector  condition  of  equation  34  yields  m  yy  scalar  equa¬ 
tions  when  expanded.  Thus,  a  solution  is  possible  as  shown  in 
Section  3.1.5  below . 


3.1.5  Representative  Iterative  Solution  For  x,^,l>  (N  given)  Via 
Steepest  Descent  First  Iteration  Tk=l) 


a . 


Given  Uj  Q,  the  nominal  trajectory,  define  the  first  iteration 

control  vector  u.  .  by 
J  >  * 


where 


A 


(35) 


(j=0,10,20,...,(N-10)) 


b. 


Then 


with 


"Fill  in"  the  control  vector  using  u.t1=u.  ,  ;i'=j ,j+l , . . .  ,j+i 

get  u.  <  where  now  i=0,l  ,2, . . .  ,N-1 . 

*  ♦  + 


j+9  to 


■+ 

x . 


i+1,1 


i=0 , 1 , . . . ,N- 1 


(36) 


-¥ 

X 


'o,l  =  X0,k  -  V  Which  is  ®iven 

Using  x^  j  solve  equation  32  for  ^  where  x^p  is  specified 
by  the  aircrew.  Entries  in  the  diagonal  matrix  [K^]  (equa¬ 
tion  25)  may  be  zero  for  control  errors  at  th end  point  that 
are  not  significant. 

Solve  equation  33  backwards  for  X-  ,  ( i=N - 1 ,  N-2 ,...,1,0) 

*  >  * 

using  the  definition  for  H.  given  by  equation  31 . 
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.1  ■  ' *  - -1 ! 


e.  Using  the  {x.  . \X.  .  ;u.  ,},  evaluate  equation  34  numerically. 

1,1  1  y  1  1 y  1 


where 


a  9H.  . 
s  =  _ JjlL 

j’1  3u .  . 

J  .1 


j  =  (0 , 10 , 20 , . . . (N- 10) ) 


Second  Iteration  (k  =  2) 

Ideally,  s.  1  should  equal  0  for  all  j.  Then  for  the  next  estimate  for 
the  optimal  control, 


where 


U  .  ,  “  X  j  s .  i 

J .  1  J .  1 


11  "  100  3H. , r 

..  J.ij 

3u . 

.)  .1 

Typically  qj  =  10%  and  J { 1  , N  j  may  be  evaluated  via  equations  22 

through  26.  The  denominator  in  the  above  is  defined  by 


3H .  2 

> 

3u .  , 

.1 . 1 


(N-10) 


3H .  , 

3u  .  , 
J  ,1 


3H  .  , 
J  > 1 

3u  .  . 
J  .  1 


Given  iu. 

'  i 


■’f ill  in"  the  control  sequence  to  obtain  the  dense  se¬ 


quence  {Uj  ,^}  for  i=0,l,2 . N-l.  Then  repeat  the  steps  of  the  first 

iteration  and  obtain  a  new  control  sequence  {u.  „}  where 

J  »  d 


22 


(40) 


Subsequent  Iterations  (k=3 .  . .  N ) 

In  general,  we  have  for  the  kth  iteration 


3H.  . 
-  J,i.k 


(41 ) 


t 


where  decreases  from  q1=10%  to  perhaps  1%  as  k  increases  according 
to  some  assumed  schedule  that  can  be  determined  by  computer  simula¬ 
tion.  Finally,  the  next  control  update  is  given  by 


(42) 


Uj,k+1  Uj,k’Xk8j,k 


(43) 


Termination  Criteria 

If  k  becomes  greater  than  k  where  k _ „  is  a  limit  that  will  be  de¬ 
max  max 

termined  by  computer  simulation,  we  terminate  the  optimal  trajectory 
generator  .lgorithm.  No  attempt  to  further  refine  N  is  then  made. 
J{k  ,N}  is,  however,  displayed  to  the  pilots  as  is  J ( 1 , N }  the  cost  of 

III  3  A 

the  nominal  trajectory.  The  pilots  may  then  decide  to  utilize  the  trajec¬ 
tory  {x.  ^  ,Uj  kmax^  a^ter  inspection  and  comparison  with  the 
nominal  trajectory, 
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where  e!  is  determined  by  the  simulation  and  where 


8H. 


U 


3u 


i,k 


( w  1 


N-l 

1 

i-0 


1 

'H 

35 

CO 

T 

[w] 

raH.  ."1 

u 

8u.  . 

L  UJ 

ai.  . 

L  UJ 

with  [w]  also  determined  via  simulation. 
Multi-Pass  Refinement  on  N 


Given  a  convergent  steepest  descent  solution  for  the  control  with 


N 


A  *0 
Ar 


repeat  the  iterative  solution  until  the  termination  criteria  is  achieved 
using  N=N2>N!  (say  10%  greater).  If  this  second  pass  is  successful, 
then  again  repeat  the  iterative  solution  procedure  using  N=N3<N!  (say 
10%  less) . 

Then  compare  the  computed  values  of  the  cost  functions  J{k1(Ni}, 

J{k2,N2}  and  J{k3,N3}to  see  which  is  least  and  choose  the  associated 

N0  and  {u.  ,  ^  }  sequence. 

*  *>•»£, 


While  an  initial  guess  at  the  control  with  N=N3  is  simply  the  control 
{u.  y}  suitably  truncated,  an  ad  hoc  approach  must  be  used  when  N=N2 
to  come  up  with  an  initial  guess  for  the  control  {u.  n}  in  the  discrete 
interval 

i  =  N ! ,  Ni+ 1 ,  N  i+2 , . . . , N 2 - 1 

A  possible  value  is  to  set  {u.  _ j  q  to  initiate  our  steepest  des¬ 

cent  iterative  procedures. 


(45) 
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After  successful  completion  of  the  complete  optimal  trajectory  determina¬ 
tion  the  {x.,u.}  sequence  is  displayed  graphically  to  the  pilots  for  their 

approval  before  being  inserted  into  the  autopilot.  The  pilots  may  elect 
to  accept  or  reject  the  trajectory  or  even  change  some  of  the  penalty 
functions  or  constraints  and  compute  a  new  optimal  trajectory. 

3.1.6  L^iution  Constraints 


The  use  of  the  penalty  function  approach  rather  than  absolute 
bounds  to  constrain  the  aircraft  state  or  to  constrain  the  control  vector 
commands  means  that  it  is  possible  that  the  original  absolute  bounds 
may  not  always  be  satisfied  within  the  optimal  trajectory.  For  crew 
safety  then,  if  the  generated  commands  exceed  the  absolute  bounds 
specified  by  the  crew  or  generated  internally  by  the  computer,  this 
condition  must  be  flagged  for  the  crew's  attention  before  the  "optimal" 
trajectory  is  engaged  into  the  flight  controller's  memory. 

The  complete  set  of  partial  derivatives  required  by  the  algorithm  will  be 
obtained  during  the  actual  algorithm  development.  Some  of  the  required 
derivatives,  however,  are  analytically  available.  For  example, 


3x 


i+1 


Arcos{i|^ } 


On  the  other  hand,  the  partial  derivatives 


3T{2i,V.,n.} 


A 


(4b) 


(47) 


must  be  evaluated  numerically  via  the  definition  of  a  partial  derivative, 
i.e. , 


A 


n.  +  6n.}  -  T{z.,Vi,n.} 

i 


using  stored  tables 

3.1.7  Trajectory  Generator  Review 

In  compliance  with  the  overall  program  plan,  System  Control, 
Inc.  (SCI)  reviewed  the  preceeding  trajectory  generator  algorithm  as 
derived  from  non-linear  optimal  control.  They  recommended  abandoning 
the  general  approach  as  described  in  favor  of  adopting  a  piecemeal  pro¬ 
cedure  for  generating  the  trajectory.  This  procedure  was 


25 


a.  To  generate  and  display  the  horizontal  projection  of  the  tra¬ 
jectory,  and 

b.  After  pilot  acceptance,  compute  the  associated  vertical  posi¬ 
tion/speed  profile  of  the  trajectory  using  assigned  altitudes 
and  minimum  fuel  algorithm. 

The  rationale  of  SCI  (Vt)  for  this  recommendation  is  stated  below: 

a.  The  general  non-linear  optimal  control  theory  approach  to 
computing  trajectories  is  too  risky  for  this  type  of  investiga¬ 
tion .  The  implementation  and  successful  development  of  this 
type  of  algorithm  will  be  time  consuming  with  a  high  possibil¬ 
ity  of  failure.  Furthermore,  the  trajectories  computed  with 
this  type  of  algorithm  will  not  be  any  more  acceptable  from  an 
operational,  flight  control  point  of  view  than  those  generated 
using  more  straightforward  techniques  In  other  words,  the 
performance  improvements  of  an  "optimal  control"  trajectory 
generator  for  the  4-D  end  point  attainment  problem  does  not 
warrant  the  development  time,  risk,  or  on-board  computer  re¬ 
sources  . 

b,  To  successfully  compute  realizable  trajectories,  the  optimal 
control  approach  may  need  to  be  modified.  The  penalty  func¬ 
tion  method  for  modeling  speed,  attitude,  and  thrust  con¬ 
straints  makes  these  constraints  "soft".  There  may  be  places 
on  the  trajectory  where  the  constraints  will  be  violated,  In 
fact,  because  of  the  characteristics  of  the  endpoint  attainment 
problem,  the  constraints  imposed  by  the  penalty  function 
method  will  probably  always  be  violated  at  the  end  of  the  tra¬ 
jectory  .  The  penalty  function  method  of  imposing  constraints 
is  the  simplest  to  implement  in  a  computer.  To  place  stiffer 
constraints  on  the  optimal  control  trajectories  significantly  in¬ 
creases  the  complexity  of  the  numerical  technique  needed  to 
compute  the  optimal  control  solution. 

In  subsequent  discussion  with  the  ,>-.ir  Force  program  officer  the  sug¬ 
gested  piecemeal  trajectory  generator  approach  was  not  implemented  be¬ 
cause  it  did  not  satisfy  the  intent  of  the  program  and  it  is  not  signifi¬ 
cantly  different  from  the  classical  trajectory  generation  technique  al¬ 
ready  developed  or.  the  IFTC  program. 

Instead,  the  program  content  was  modified  to  use  the  existing  classical 
trajectory  generation  technique  for  the  simulation,  but  increasing  the 
capabilities  of  these  algorithms  for  the  survivability  of  military  trans¬ 
port,  bomber  vehicles  in  the  presence  of  SAM  and  AAA  site  ground 
threats . 
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The  general  task,  as  described  in  Section  3.3,  was  assigned  to  SCI  as 
part  of  their  analytic  support  contract.  SCI  was  qualified  for  this  task 
because  of  their  previous  development  :>f  minimum  exposure  profile  gen¬ 
erators  on  the  Navy-funded  Advanced  Weapons  Management  System  pro¬ 
gram  . 

3.2  OPTIMAL  CONTROLLER  ALGORITHM  DEVELOPMENT 

The  general  procedure  of  developing  the  linear  quadratic  optimal 
control  algorithm  included  the  following  steps: 

a.  Simplification  of  the  aerodynamic  coefficients  for  the  KC-135 
aircraft. 

b .  Generation  and  simplification  of  the  aircraft  state  equations 
using  the  developed  coefficients  of  step  a. 

c.  Linearization  of  the  state  equations  in  order  to  apply  the 
linear  quadratic  estimator. 

d.  Design  of  the  linear  quadratic  estimator  by  selection  of  the 
state  and  control  vectors  of  the  cost  functional  and  analysis 
of  the  system  response  to  deterministic  inputs  and  sensitivity 
to  aircraft  flight  dynamics. 

e .  Investigation  of  mechanization  requirements . 


3.2.1  Simplified  Aerodynamic  Coefficients 


A  primary  objective  in  developing  the  state  model  for  the  (LQ) 
regulator  is  to  obtain  the  simplest  aircraft  model  that  still  retains  an 
adequate  description  of  the  airplane  dynamics.  In  order  to  achieve  this 
objective,  the  following  assumptions  were  made: 


□  The  aircraft  is  rigid  body  so  that  aeroelastio  effects  were  not 
modeled. 


□  The  aircraft  mass,  moment  and  products  of  inertia,  and  cen¬ 
ter  of  gravity  are  constant. 

□  Ground  effect  and  landing  gear  coefficients  were  not  modeled 
since  landings  were  not  simulated. 

□  Refuel  boom  coefficients  are  negligible. 


a  Spoilers  and  flaps  are  not  deployed. 
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□  Stabilizer  angle  is  fixed  at  -6°. 

□  Ice  doesn't  form  on  wings. 

In  addition,  engineering  judgement  was  used  to  simplify  the  remaining 
significant  dimensionless  aerodynamic  force  and  moment  coefficients.  In 
most  cases  the  resulting  coefficient  is  within  an  order  of  magnitude  of 
the  actual  value. 


Appendix  A  of  this  report  contains  the  KC-135  aerodynamic  force  and 
moment  coefficients  as  they  appear  in  Boeing  Document  No.  D3-9781-2. 
Many  of  the  components  of  the  force  and  moment  coefficients  are  elimi¬ 
nated  as  a  consequence  of  the  preceding  simplifying  assumptions. 

The  simplified  dimensionless  aerodynamic  coefficients  (Cn,  Cv>  CT  ,  C0, 
Cm»  Cn)  are  then  defined  as, 


where  CT  ,  CU,  C  ,  C0,  C  ,  and  Cv  are  the  aerodynamic  coefficients 
ij  in  z  n  x 

of  list,  drag,  pitch-moment,  roll-moment,  yaw-moment,  and  side  forces, 
respectively.  For  ihese  coefficient  equations,  the  term  definitions  are: 


p  =  airplane  sideslip  angle  (deg) 

6^  =  rudder  deflection  angle  (deg) 

P  =  airplane  sideslip  angle  rate  (rad/sec) 

a  =  airplane  angle  of  attack  relative  to  the 

wing  chord  plane  (deg) 


p 

roll  rate 
(rad/sec) 

about 

the 

stability 

axis 

r 

= 

yaw  rate 
(rad/sec) 

about 

the 

stability 

axis 

6„ 

- 

aileron  deflection  angle  (deg) 

q 

- 

pitch  rate 
(rad/sec) 

about 

the 

stability 

axis 

a  =  airplane  angle-of-attack  rate  (rad/sec) 


6e  =  elevator  deflection  angle  (deg) 

M  =  Mach  number 

b  =  wing  span  =  130.83  feet  for  the  KC-135 
aircraft 

V  =  true  airspeed  (ft/sec) 

c  =  wing  mnan  aerodynamic  chord  =  20.16 

feet  for  the  KC-135  aircraft 

Additional  simplifications  to  these  coefficients  can  be  made  by  consider¬ 
ing  the  dynamics  of  the  lateral  and  longitudinal  modes  of  an  airframe.3 
The  approximate  damping  factor  and  natural  frequency  for  the  quadratic 
modes  and  time  constants  for  the  first  order  modes  are 


(CT  +  C,J  -  2u  C 
v  L  Ir  cm 


3Dynamics  of  Airplane,  Northrup  Corporation,  Aircrafts  Division , 
Document  number  AE-61-4-11. 


short 

period 


I  / _i_\ 
T  V'yy) 
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,k 

yy  V 


Examination  of  equations  49  through  56  reveals  that  the  following  coef¬ 
ficients  do  not  appear. 

CT  ,  CT  ,  C  ,  and  C 
L*  ’  L  ni’  n 

a  q  P  P 

These  coefficients  were  eliminated  because  they  were  not  deemed  signifi¬ 
cant  to  longitudinal  or  lateral  mode 
characteristics.  The  effects  of 

C,  ,  C„  ,  and  C, 


are  small  (£  0. 004/degree)  and  were  also  eliminated.  If  the  lesulting 
dimensionless  aero  force  and  moment  equations  are  used  in  the  state 
mode,  the  term 


•i  (°0 


in  the4  dimensionless  pitching  momemt  equation  couples  the  state  deriva¬ 
tives  a  and  q  .  This  term  cannot  be  neglected  as  it  contributes  rough¬ 
ly  20%  of  thes  short  period  mode  damping.  It  can  be  shown,  however 
that  modifying  the  term 


W 


to  (C+C  )  §£. 

ma  Kq  2VT 


and  deleting  the  term 


■i("  i) 


will  decouple  the  state  derivatives  and  not  alter  the  short  period  mode 
damping.  These  further  simplifications,  with  substituted  into  the 
Cq  equation,  reduce  the  dimensionless  force  and  moment  coefficients  of 
equation  48  to, 
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From  the  graphs  in  Appendix  A,  the  simplified  values  of  terms  in  these 
coefficients  for  the  KC-135  aircraft  are: 


.0796  (1/deg) 


.0191  (1/deg) 


-5.  (l/rad) 


3.2.2  Aircraft  State  Model 
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A  complete  six-degree-of-freedom  aircraft  state  model  was  de¬ 
rived,  as  defined  in  Section  3.4,  Simulator  Development,  for  closed- loop 
simulation  of  the  optimal  control  algorithms  using  the  transport  cockpit. 

These  equations  are, 
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For  equations  58  and  59  the  term  definitions  are 


P 

a 


VVfs 


PS’Vrs 


p 

a 

<M.4> 

XW,YW’ZW 


X  , Y  ,Z, 
s  s  s 


Ws 


T 


W 

M 


true  airspeed  rate  (ft/sec2) 

sideslip  angle  rate  (rad/s ec) 

angle-of-attack  rate  (rad/sec) 

roll,  pitch,  and  yaw  accelerations  about 
stability  axis  (rad/sec*) 

roll,  pitch,  and  yaw  rates  about 
stability  axis  (rad/sec) 

true  airspeed  (ft/sec) 

sideslip  angle  (rad) 

angle  of  attack  (rad) 

roll,  pitch,  yaw  attitude  rate  (rad/sec) 

roll,  pitch,  yaw  attitude  (rad) 

external  forces  along  the  wind  axes 
(lbs) 

aerodynamic  forces  along  stability  axes 
(ibs) 

aerodynamic  moments  about  stability 
axes  (ft-lbs) 

engine  thrust  (lbs) 

aircraft  weight  (132,146  lbs  for  KC- 135 
aircraft) 

aircraft  mass  (4104  slugs  for  KC-135 
aircraft) 


U  >  h » h » Is  >  1 10 


inertial  constants  defined  in  Section 
3. 4. 1.2  of  the  roll,  pitch,  yaw  and 
product  moments  of  inertial  of  the 
aircraft 
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S  =  wing  reference  area  (ft2) 

c  =  mean  aerodynamic  chrod  (ft) 

b  =  wing  span  (ft) 

M  =  Mach  number 

p  =  air  density  (slugs/ft3) 

Cn.Cv.CT  ,C0,C  ,C  =  coefficients  of  aircraft  aerodynamics  as 

defined  in  Section  3.2.1 

By  substitution  of  the  aerodynamic  coefficients  of  Section  3.2.1  and  the 
above  definition  equations  and  noting  that  the  angle  of  attack  with 
respect  to  the  wing  chord  plane  is  related  to  the  angle  of  attack  of  the 
x  body  axis  by  the  equation  awrp  =  a  +  2°,  the  final  state  model  equa¬ 
tions  are, 
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To  reduce  the  complexity  of  the  required  solution  for  the  linear  quadra¬ 
tic  optimal  controller,  these  six-degree-of-freedom  equations  were  re¬ 
duced  to  simple  point  mass  equations  of  motion. 

The  vector  equation  describing  the  motion  of  the  center  of  mass  of  an 
airplane  with  respect  to  an  inertial  frame  is 


F  = 

d  ,, 
m  tt  V 

dt  - 

where 

F 

external  force  acting  upon 

the 

airplane 

m  = 

mass  of  airplane 

V 

velocity  of  the  center  of 

mass 

of  the 

airplane  with  respect  to 

an 

inertial 

frame 

The  inertial  reference  frame  is  fixed  to  the  earth.  The  z  axis  is  along 
the  gravity  vector.  The  x  axis  points  true  north  and  the  y  axis  com¬ 
plete  a  right-handed  orthogonal  coordinate  system.  The  total  velocity 
vector  (V)  of  the  aircrr.:  with  respect  to  ground  is  equal  to  the  vector 
sum  of  the  aircraft  velocity  with  respect  to  air  (Yair^  and  the  w^d 

velocity  (V^)  the  velocity  of  the  air  mass  with  respect  to  inertial  co¬ 
ordinates,  Rewriting  equation  61 
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m  (dt  -AIR  +  dt  V 


(62) 


The  equations  of  motion  are  written  with  respect  to  an  Eulerian  axis  sys 
tern.  This  coordinate  system,  called  the  wind  axis  system,  is  fixed  to 
the  center  of  mass  of  the  airplane.  The  x  axis  is  coincident  with  the 
total  velocity  of  the  airplane.  The  z  axis  lies  in  the  aircraft  plane  of 
symmetry  perpendicular  to  the  x  axis  and  positive  downward.  The  y 
axis  completes  a  right-handed  orthogonal  system.  If  W  is  the  rotational 
rate  of  the  wind  with  respect  to  inertial  coordinates ,  then  by  rules  for 
differentiation  of  vector, 
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and  Rj^  is  the  direction  cosine  matrix  relating  inertial  to  wind  axis  co¬ 
ordinates  .  The  external  forces  acting  on  the  airplane  consist  of  grav¬ 
ity,  aerodynamic  forces  and  thrust.  These  forces  resolved  along  the 
wind  axis  system  are 
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The  term  definitions  of  equation  64  are 

D  =  drag 

g  =  acceleration  of  gravity  constant 

T  -  thrust 

a  =  angle  between  thrust  vector  and  velo 

city  vector 

Y  =  flight  path  angle 


4>  =  bank  angle 

L  =  lift 


Considering  only  the  horizontal  component  of  wind  in  the  ground  frame 


[VyCos  4<w  Vwsin4«w  0) 


(07) 


where  Vw  is  the  magnitude  of  the  horizontal  wind  and  4\y  is  the  wind 
direction  in  inertial  coordinates.  The  direction  cosine  matrix 
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is  used  to  transform  the  time  derivative  of  wind  vector  ^  from  the 

inertial  reference  frame  to  the  wind  axis  frame.  Substituting  equations 
64  through  68  into  equation  63, 


Lift  is  modeled  by  the  relation 

L  =  3sp(h)VzSCL 

and  drag  by  the  equation 
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The  dimensionless  aerodynamic  coefficients  C.  and  Cp.  were  derived  in 
Section  3.2.1  as 
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where 


a  +  2‘ 


The  air  density  p(h)  is  approximated  by  the  quadratic 


P0 [ 1  ‘  '  1,835(65536)  +  1 


where  is  the  air  density  at  sea  level. 

The  aircraft  position  relative  to  earth  is  obtained  by  resolving  the  air¬ 
craft  velocity  vector  along  the  inertial  reference  frame  axis. 


(7',) 


To  steer  the  aircraft  in  the  horizontal  plane,  it  is  desirable  to  have 
along-track  and  cross-track  distances  available  as  state  variables.  Re¬ 
ferring  to  Figure  2,  a  coordinate  frame  is  defined  with  its  origin  at  a 
reference  position  on  the  nominal  horizontal  trajectory  such  that  the  dis¬ 
tance  between  the  reference  and  aircraft  position  is  minimized.  De¬ 
fining  the  reference  position  in  this  manner  ensures  that  the  control 
corrections  will  always  be  steering  the  aircraft  in  order  to  maintain  the 
smallest  spatial,  errors  possible.  Along-track  and  cross-track  rates  are 
obtained  from  x  and  y  using  the  relation 


(76) 
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Modifying  equation  74 
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MOVING  TARGET  COORDINATE  SYSTEM 
FIGURE  2 


The  horizontal  wind  magnitude  (Vw)  and  wind  direction  (t^)  will  be 

modeled  as  deterministic  disturbances.  The  disturbance  dynamics  are  in 
the  form  of 


where 


where  x  consists  of  a  nominal  state  vector  (Xq)  and  a  state  perturba¬ 
tion  vector  (fix)  and  U  consists  of  a  nominal  control  vector  (U^)  and  a 
control  correction  vector  (6U). 
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where,  in  the  general  form 
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The  partial  derivatives  that  comprise  the  non-zero  elements  of  the  ma¬ 
trices  A  and  B  are  obtained  by  performing  the  indicated  operation  on 
equation  79  and  evaluating  at  x=xq,  u=u , .  Then, 
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Nominal  horizontal  wind  magnitude  and  direction  gradients  were  assumed 
to  be  equal  to  zero. 

3.2.4  Design  of  Linear  Quadratic  Regulator 


The  state  model  developed  in  this  previous  section  may  be  writ¬ 
ten  in  the  form, 
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6a  =  difference  between  measured  and  nominal 
angle-of-attack  (rad) 

6<|>  =  difference  between  measured  and  nominal 

bank  angle  (rad) 

6T  =  difference  between  measured  nominal 
thrust  (lbs) 
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H(t) 

F(t) 

B(t) 


6x6  plant  system  matrix 
6x2  disturbance  input  matrix 
2x2  disturbance  system  matrix 
6x3  plant  input  matrix 


The  form  of  these  matrices  is 
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and  the  explicit  form  of  the  non-zero  elements  are  the  partial  deriva¬ 
tives  of  the  state  equations,  evaluated  at  x_Uo,  as  shown  in  the  pre¬ 
vious  section. 
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The  quadratic  cost  functional  to  be  minimized  is 


06  t  t  few  oirxctn  T 

J(t)  =  J  {[x(t)1  z(t)J]  j^0  oJj_zCt)J+  u1  (t)R(t  )u(t ) }  dt  (89) 

Note  that,  since  the  control  correction  vector  cannot  influence  distur¬ 
bances  z(t),  these  disturbances  are  not  included  in  the  cost  function. 

This  steady-state  linear  quadratic  problem  has  a  solution  provided  (1) 

[A,B]  is  controllable,  and  (2)  F  is  strictly  stable  (eigenvalues  in  left 
half  plane). 

The  optimal  control  law  consists  of  plant  state  feedbacks  and  distur¬ 
bance  state  feedbacks,  i.e., 

U(t)  =  -GX(t)  -  PZ(t)  (90) 

This  is  shown  in  block  diagram  form  in  Figure  3. 


LINEAR  QUADRATIC  CONTROLLER  BLOCK  DIAGRAM 

FIGURE  3 
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In  this  diagram, 

G  is  a  3x6  constant  plant  state  gain  matrix  and  is  given  by 

G  =  R-1BTK  (91) 

where  K  is  6x6  and  is  the  solution  to  the  algebraic  Ricatti  equation 

0  =  -KA  -  ATK  -  Q  +  KBR~1BTK  (92) 

Observe  that  the  gain  matrix  G  is  independent  of  the  disturbance  dyna¬ 
mics  . 

P  is  a  3x2  constant  disturbance  state  gain  matrix  and  is  given  by 

P  =  R_1BTM  (93) 

where  M  is  6x2  and  is  the  solution  to 

0  =  -KH  -  ATM  -  MF  +  KBR"1BTM  (94) 

The  solution  of  these  equations  is  dependent  on  the  values  of  Q(t)  and 
R(5).  Assume  that 

Q(t)  -  diag[Qu(t)  Q2 2 ( t)  Q3 3 (t)  Q44 (t)  Q5g(t)  Qe6(5)l 

(95) 

R(t)  =  diag[R^j(t)  Rgj(t)  R33(t)] 


The  weighting  matrices  were  chosen  to  be  diagonal  since  this  allows  the 
states  and  controls  to  be  penalized  individually. 

Selection  of  the  weighting  matrices  §  and  R  was  accomplished  in  four 
steps . 

Step  1:  Obtain  a  good  first  guess  of  the  weighting  matrices  § 

and  R  for  a  specific  flight  condition. 

Step  2:  Assuming  the  same  flight  condition  as  in  step  1,  indi¬ 

vidually  vary  the  components  of  the  weighting  matri¬ 
ces  to  determine  the  effect  on  the  resulting  closed 
loop  system  eigenvalues. 

Step  3:  Using  the  weighting  matrices  selected  in  step  1, 

study  the  locus  of  the  closed  loop  system  eigenvalues 
as  a  function  of  flight  condition. 
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Step  4: 


Use  the  results  of  steps  2  and  3  to  schedule  the 
weighting  matrices  §  and  R  such  that  critically 
damped  eigenvalues  are  obtained  for  varying  flight 
conditions . 

3. 2.4.1  LQ  Design  -  Step  1 

A  common  method  of  initially  selecting  the  Q  and  R  matrices  is 
to  assign  maximum  values  to  the  plant  states  and  maximum  allowable  con¬ 
trols  corresponding  to  those  values.  This  information  is  communicated 
to  the  mathematics  by  restructuring  the  cost  functional  in  the  form 


J  = 


where 


(96) 


The  maximum  values  assigned  to  the  state  and  control  perturbations  are 
shown  in  Table  I  for  a  wings  level  aircraft  flight  condition  having  a  low 
dynamic  pressure  q  corresponding  to  an  altitude  of  5000  feet  and  a  velo¬ 
city  of  274.8  feet  per  second, 


The  structure  of  the  optimal  feedback  control  corresponding  to  the  solu¬ 
tion  of  several  different  steady-state  linear  quadratic  problem  formula¬ 
tions  is  given  in  Table  II.  For  wings  level  and  no  wind  disturbances 
the  lateral  control  6<Kt)  is  independent  of  the  vertical  (6y(t),  6h(t)) 
and  longitudinal  (6Vlt),  6x^^(t))  .node  state  feedbacks.  The  lateral 

mode  is,  therefore,  decoupled  from  the  longitudinal  and  vertical  modes. 
Coupling  exists,  however,  between  the  longitudinal  and  vertical  modes. 
For  a  nominal  aircraft  bank  angle  of  45°,  all  three  path  modes  are 
coupled.  When  wind  disturbances  are  included,  the  disturbance  states 
begin  to  appear  in  the  control  equations.  For  a  50  ft/sec  tailwind,  the 
wind  magnitude  perturbation  contributes  to  the  perturbed  angle-of- 
attack  and  thrust  commands,  The  wind  direction  perturbation  contri¬ 
butes  tc  the  perturbed  sank  angle  command.  For  a  50  ft/sec  cross- 
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TABLE  1 

INITIAL  g  AN1)  H  SELECTION 


6V 


max 


641 


max 


6V, 


max 


50  ft/sec 

7.5° 

5° 


6xat 

=  2000 

ft 

6a 

max 

II 

C n 
0 

max 

6yCT 

=  1000 

ft 

66 

Tmax 

=  7.5° 

max 

6h 

max 

100  ft 

6T 

max 

=  12000  lbs 

TABLE  II 

SOLUTION  OF  STEADY  STATE  LQ  PROBLEM 


•  NO  WIND  DISTURBANCES  -  LOW  q  FLIGHT  CONDITION  WITH  NOMINAL  PERFORMANCE  INDEX 

•  NINOS  LF.Vl-l. 

A»U)  •  MH)409AV(t)  4  1.236Y(I>  -  .  0000 l  3A  xAT  ( t  )  I  .00087/Ah(t) 

.'4(0  -  1.27A*(t>  *  .0001ilAyCT(t) 

AT  ( t )  -  379.6VU)  4  15302.  Ar(t)  4  5.?16xA.f(0  I  36.4Ah(0 

•  43°  ROLL  ATTITUDE 

An  (0  -  0003176V(t)  4  , 31244(0  +  l.G436y(0  -  .00001424xAT(O 

4-  .0000368AyCT(t)  +  .000?326h(t) 

64(0  -  -0.00129AV(t)  +  1.05244'(O  -  1 . 4B4y (O  4  . 0000085ixAT(r ) 

4  .0001194yCT(O  -  .0II052(i6h(t) 

AT  <  t  >  -  41 . 8AV(t)  4  316 ,44*Ct)  4  28346,  4rU)  4  5.62AxAT(0 
4  ,2524yCT<t)  +  41,6Ah<0 

•  WITH  WIND  DISTURBANCES  -  LOW  5  FLiCHT  CONDITION  WITH  NOMINAL  PERFORMANCE 
INDEX  AND  WINGS  LEVEL 

•  30  F r/SEC  TAILWIND 

Au (l )  -  ,000409AV(t)  4  l.236y(t)  -  . 00001 36xAT(t)  4  ,0008276h(t) 

-  .  0002946Vy(t) 

44U)  *  1.2764(c)  4  ,0001316yCT(t)  4  ,173i*u(0 

6T(t)  -  379  .  AV(t)  4  15302. 6y(c)  4  5 . 71AxA„.(0  4  36.46h(t) 

4  304 . 5Vy(t ) 

•  50  FT/ SEC  CRGSSWIND 

6a ( t )  -  . 0004096V(t)  4  i.234y(t)  -  . 0000136xAT(t)  4  ,0008276h(t) 

-  . 0000000 346Vy(t)  4  ,  008464u(O 

A*U)  •  l .  2764 (t)  y  ,00013l4yCT(t)  4  ,002236Vu(t)  4  ,0000264u(O 

6T(  t )  -  379 . 6V(  t )  4  15302. 6y(c)  4  5 . 71  Ax^U )  4  3t  46tii  ) 

4  . 035i4Vy(C)  -  12400. 64w(t) 


wind,  the  disturbance  state  feedback  gains  change  so  that  the  wind 
magnitude  perturbation  affects  mainly  the  perturbed  bank  angle  com¬ 
mand  while  the  wind  direction  perturbaLicn  affects  mainly  the  perturbed 
angle-of-attack  and  thrust  commands 

3. 2.4. 2  LQ  Design  -  Step  2 

An  eigenvalue  analysis  was  made  to  determine  the  effects  vary¬ 
ing  individual  components  of  the  Q  and  R  matrices  had  on  the  stability 
of  the  optimal  controller.  Figures  4  through  12  are  root  locus  type  dia¬ 
grams  which  illustrate  these  effects.  The  following  conclusions  can  be 
drawn  from  these  figures. 

□  Three  sets  of  eigenvalues  can  be  identified  as  corresponding 
to  the  longitudinal,  lateral,  and  vertical  path  modes, 

□  Longitudinal  path  mode  has  a  bandwidth  of  the  order  of 
.035  rad/sec. 

□  Lateral  path  mode  has  a  bandwidth  of  the  order  of  ,06 
rad/sec. 

□  Vertical  path  mode  has  a  bandwidth  of  the  order  of  .5 
rad/sec . 

□  Decreasing  decreases  the  longitudinal  mode  bandwidth 

and  also  effects  the  vertical  mode  bandwidth. 


□ 

□ 

□ 

□ 

n 

u 


Decreasing  6t|i  decreases  the  lateral  mode  bandwidth, 
max 

Decreasing  6y„  ,  decreases  the  vertical  mode  bandwidth . 
°  ‘max 

Increasing  6x^,p  decreases  the  longitudinal  mode  band- 
widith.  1  max 


Increasing 


6yCT 


max 


decreases  the  lateral  mode  bandwidth . 


Increasing  6h  decreases  the  vertical  mode  bandwidth, 
max 

Increasing  6a  ^  decreases  the  vertical  mode  bandwidth  and 
also  effects  the  longitudinal  mode  bandwidth. 


□  Increasing  6d>  decreases  the  lateral  mode  bandwidth. 

max 

□  Increasing  6T  decreases  the  longitudinal  mode  bandwidth. 

max 
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3. 2. 4. 3  LQ  Design  -  Step  3 

A  sensitivity  analysis  was  performed  to  study  the  variation  of 
the  closed-loop  eigenvalues  as  the  nominal  flight  path  angle  (Yo)>  nomi¬ 
nal  bank  angle  ( 4>0 ) ,  nominal  altitude  (h0),  or  nominal  velocity  (V0) 
varies.  Nominal  thrust  (T0)  and  angle-of-attack  (a0)  were  adjusted  so 
that  the  aircraft  is  trimmed,  Root  locus  type  pilots  illustrating  the  re¬ 
sults  of  this  analysis  are  shown  in  Figures  13  through  16,  The  follow¬ 
ing  observations  can  be  made: 

n  Varying  the  nominal  flight  path  angle  has  negligible  effect  on 
path  mode  dynamics . 

n  The  vertical  path  mode  dynamics  are  strongly  influenced  by 
varying  nominal  altitude  and  to  a  lesser  extent  by  varying 
nominal  velocity . 

□  The  lateral  path  mode  dynamics  are  strongly  affected  by  vary¬ 
ing  nominal  velocity . 

□  The  longitudinal  path  mode  dynamics  are  only  slightly  affected 
by  varying  nominal  altitude,  velocity,  and  bank  angle, 
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ROOT  LOCUS  TOR  NOMINAL  ALTITUDE  VARIATION 


.  NOMINAL  FUOHT  CONDITION  WITH  ALTITUDE  CHANGES 


ROOT  LOCUS  FOR  NOMINAL  VELOCITY  VARIATION 

FIGURE  16 
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3. 2.4.4  LQ  Design  -  Step  4 


The  weighting  matrices  Q  and  R  are  scheduled  as  a  function  of 
flight  condition  by  selecting  a  performance  index  for  a  sample  of  re¬ 
presentative  flight  conditions  that  yielded  eigenvalues  that  are  approxi¬ 
mately  critically  damped.  The  rationale  for  critically  damped  eigen¬ 
values  is  to  minimize  overshoots  which  are  undesirable  in  path  guid¬ 
ance.  The  performance  index  was  modified  by  changing  the  maximum 

velocity  error  (6V  ),  heading  error  (6ip  ),  and  flight  path  angle 

HiaX  max 

error  (6y  )  terms.  Figures  17,  18,  and  19  show  the  error  schedul- 

max 

ing  as  a  function  of  velocity,  altitude,  and  bank  angle.  Table  Ill  pre¬ 
sents  the  final  weighting  matrix  selection  illustrating  the  actual  imple¬ 
mentation  of  the  6V  ,  6di  ,  and  6v  scheduling. 

max  max  'max 

3. 2.4.5  Mechanization  Requirements 

The  design  of  the  optimal  controller,  as  defined,  incorporates 
two  principal  approximations  to  keep  the  requirement  for  airborne  pro¬ 
cessing  capability  within  reason.  These  restrictions  were  (a)  the  use 
of  a  point  mass  aircraft  model,  and  (b)  the  implementation  of  a  propor¬ 
tional  control  system  that  excluded  integral  terms.  These  simplified  mo¬ 
dels  require  the  scheduling  of  the  weighting  matrices,  Q  and  R,  as  a 
function  of  aircraft  flight  condition  (as  previously  shown)  and  analysis 
of  the  inner  loop  gains  of  the  normal  Stability  Augmentation  System 
(SAS)  of  the  KC-135  aircraft. 

In  addition,  simplified  aircraft  control  systems  required  either  angle-of- 
attack  or  normal  acceleration  control  signals  for  the  outer  loop  pitch 
axis.  This  type  of  control  system  is  considerably  more  difficult  to  im¬ 
plement  than  a  conventional  pitch  attitude  control  system.  In  fact, 
without  the  integral  terms  in  the  control  system,  the  angle  of  attack  or 
normal  acceleration  control  signals  in  the  outer  loop  would  also  have  to 
be  scheduled  in  gain,  as  a  function  of  flight  mode,  to  maintain  desired 
stability . 

Performance  improvements  of  an  optimal  control  system  compared  to  a 
conventional  control  system  on  a  KC-135  aircraft  would  be  difficult  to 
measure  if  :.he  only  criteria  were  the  arrival  at  a  rendezvous  point  or 
CARP  at  a  given  time.  However,  some  of  the  negative  comparisons  may 
be  more  apparent  to  the  air  crew.  The  crossfeed  terms  derived  by  the 
optimal  control  system  would  undoubtedly  provide  smaller  profile  errors, 
but  the  transient  responses  to  pilot  overrides  and  at  profile  boundaries 
when  new  gains  were  generated  may  be  very  disconcerting  t.c  pilots. 
Also,  the  crossfeed  terms  to  the  throttle  would  probably  cause  an  exces¬ 
sive  amount  of  throttle  activity  . 
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These  mechanization  problems  resulted  in  a  decision  by  the  AFKDL  pro¬ 
gram  manager  to  terminate  the  optimal  controller  development  task  and 
proceed  with  the  effort  of  demonstrating  the  improved  threat  avoidance 
conventional  algorithms  and  specifying  a  preferred  method  to  demon¬ 
strate  the  integrated  flight  trajectory  system  in  a  future  flight  test 
demonstration  program. 

3.3  THREAT  AVOIDANCE  ALGORITHM 

As  defined  in  Section  3.1.7,  the  task  of  defining  an  optimal  trajec¬ 
tory  generator,  because  of  the  complexity  involved,  was  deferred  in  fa¬ 
vor  of  upgrading  the  existing  conventional  algorithm  trajectory  genera¬ 
tor  to  include  threat  avoidance  capability.  Specifically,  the  function  of 
the  algorithms  is  to  compute  mission  trajectory  modifications  for  avoiding 
or  minimizing  the  aircraft  exposure  to  ground  threats  (SAMs  or  AAA) 
during  tactical  airlift  or  transport/bomber  missions.  A  secondary  objec¬ 
tive  is  to  develop  threat  evasion  algorithms  which  compute  guidance/ 
control  cues  or  commands  for  evasive  maneuvers  after  a  threat  has  been 
launched,  detected,  and  identified. 

These  algorithms  will  be  simulated  and  evaluated  with  a  cockpit,  and 
hybrid  computer  simulator  that  includes  a  KC-135  cockpit  and  aerodyna¬ 
mic  model  plus  a  tactical  flight  management  system.  This  flight  manage¬ 
ment  system  consists  of  a  tactical  situation  display,  an  alphanumeric 
display,  and  keyboard  installed  in  the  cockpit,  and  a  4-D  trajectory 
generator  for  path  guidance  implemented  in  the  hybrid's  digital  compu¬ 
ter.  In  the  simulation,  the  increased  survivability  algorithms  will  be 
integrated  with  the  existing  4-D  trajectory  generation  algorithms,  which 
have  been  developed  on  the  USAF  Integrated  Flight  Trajectory  Control 
(IFTC)  program. 

3.3.1  Tas k  D e f ini t io n 

The  specific  scope  of  the  revised  task  included  in  analytical  char¬ 
acterization  of  the  ground-located  threats  and  definition  of  increased 
survivability  trajectory  generator  algorithms. 

As  a  background  for  the  algorithm  definition  task,  a  realistic  ground 
threat  environment  foi  tactical  airlift  or  bomber  missions  in  the  1980- 
1990  European  threatre  was  requested.  The  ground  threat  would  then 
be  analytically  modeled  with  equations  that  could  be  used  to  adequately 
represent  the  threats  on  a  tactical  situation,  display.  The  threat  avoid¬ 
ance,  minimum  exposure,  and  threat  evasion  algorithms  that  are  devel¬ 
oped  would  be  compatible  with  the  IFTC  trajectory  generator  equations 
and  the  flight  control  system  capabilities  of  the  simulated  KC-135  air¬ 
craft. 
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The  final  engineering  report  on  this  task,  as  generated  by  SCI  (Vt),  is 
contained  in  Sections  3.3.2  through  3.3.4.  This  report  briefly  reviews 
the  concept  of  resource  management  for  aircraft  survival,  threat  model¬ 
ing,  and  threat  avoidance  algorithms. 


Resource  Management  Concepts 


The  problem  of  threat  avoidance  is  best  considered  in  the  more 
comprehensive  framework  of  resource  management.  By  normalizing  the 
current  military  value  of  an  aircraft  in  terms  of  aircraft  (e.g.,  2.7  air¬ 
craft)  as  a  function  of  the  vector  of  allocable  resources,  R,  for  multiple 
objectives ,  it  can  be  shown  that 


“N-l  mN+l 

+  ...  +  Tm  II  S.tL  +  n  S.  (98) 

» j=1  ^  •> 

where 

T.  =  value  of  the  ith  objective  (i  =  1,2, 

1  ....  N) 


M.  =  probability  of  accomplishing  ith  objec- 
1  '  tive  conditioned  upon  arriving  at  ith 

objective 


m. 

l 

index  of 
objective 

interval  incorporating  ith 

S. 

J 

= 

probability 

of  surviving  jth  interval 

Vl 

(11  S. 

= 

the  probability  of  survival  to  the 

j=l  J 

aircraft  ra 

turn  base) 

The  primary  resources  available  to  an  aircraft  fall  into  three  categories. 
The  aircraft  maneuvers,  specified  by  the  three-dimensional  position  and 
velocity  as  a  function  of  time,  comprise  the  first.  The  second  is  non¬ 
expendable  resources ,  which  can  be  allocated  at  each  time  of  mission 
without  "using  up"  the  resource.  For  example,  ECM  RF  power,  fire 
control,  and  surveillance  sensors,  radar  cross  section  presented  to 
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threats  (which  ...  function  of  relative  a„ui.ude),  and  computer  capac¬ 
ity.  The  third  resource  category  includes  expendable  resources  which, 
once  used,  are  no  longer  available.  These  include  flares,  decoys, 
guns,  defensive/offensive  missiles,  etc. 

The  resource  management  appi  .  block  diagrammed  in  Figure  20,  is 
to  formulate  the  problem  in  recu  ,ive  form  and  maximize  the  index  for¬ 
mulated  in  equation  98  over  the  allowable  flight  control  subject  to  fuel 
constraints . 

R-dar  cross  section  is  controlled  by  .  '"'sk.oring  the  effect  of  aspect 
angle  and,  hence,  cross  section,  on  tut  probability  of  detection  by  var¬ 
ious  radars  along  the  flight  path,  These  probabilities  are  then  mapped 
into  the  probability  of  aircraft  survival  at  each  grid  point  in  a  three- 
dimensional  position /velocity  space.  A  similar  approach  is  used  to  model 
the  effects  of  terrain  masking.  Part-;  of  the  space  that  are  masked 
from  threat  radars  by  terrain  are  assigned  a  unity  probability  of  sur¬ 
vival,  while  portions  of  space  below  ground  level  are  given  a  probability 
of  survival  of  zero. 

As  car  be  envisioned,  the  objective  is  to  maximize  the  probability  of 
surv.val  along  all  possible  flight  palls  P-  and  returning  from  the  objec¬ 
tives,  match  the  resulting  solutions  to  boundary  conditions,  and  then 
choose  the  one  flight  path  giving  the  highest  current  military  value  of 
the  aircraft  using  models  that  describe  the  effect  on  the  performance 
measure  of  allocating  a  particular  resource  to  a  given  threat .  An  impor¬ 
tant  issue  in  the  optimal  allocation  of  resources  if  the  time-line  trade-off 
between  an  allocation  of  resources  at  a  current  time  during  a  mission 
and  the  allocation  of  the  remaining  resources  at  later  times  during  that 
mission,  The  algorithm  must  include  the  capability  for  real-time  alloca¬ 
tion  of  resources  as  threats  actually  materalize  (e.g.,  10  seconds)  while 
still  r  serving  the  necessary  resources  for  allocation  in  later  mission 
phases.  The  allocation  for  later  phases  of  doe  mission  usually  can  be 
carried  ^ut  on  a  slower  time  scale.  These  two  "fast"  a,..i  "slow"  a’loca- 
t.ion  regimes  are  not  independent,  and  a  consistent  global  optimization 
mi  be  performed;  the  expendables  are  allocated  globally,  while  the 
non  expendables  may  be  allocated  locally. 

One  possible  approach  to  the  problem  taken  by  SCI4  is  a  modificition  of 
the  dynamic  programming  prccedure  for  solving  multistage  decision  pro¬ 
cesses  to  make  the  procedure  sur  lble  for  airborne  computation. 

4,).P.  March  and  M.  Grossberg,  "Advanced  Weapons  Management 
Systems",  Pa:  Tic  Missile  Test  Center,  April  1978. 
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3.3.3  Threat  Avoidance  Review 


The  implementation  of  threat-avoidance  algorithms  can  be  per¬ 
formed  at  various  levels  of  sophistication .  The  algorithm  selected  is 
directly  related  to  the  quality  of  the  threat  models . 

Using  a  cylindrical  threat  volume  around  a  known  threat  size,  algo¬ 
rithms  developed  by  LSI  on  the  IFTC  simulation  skirt  around  the  peri¬ 
phery  of  the  threat  volume  using  a  trajectory-generation  technique  link¬ 
ing  a  set  of  dynamically  generated  waypoints.  However,  the  approach 
is  not  satisfactory  in  multi-threat  environments  where  the  threat  vol¬ 
umes  are  overlapping  and  the  deployment  of  expendable  resources  is  not 
accounted  for. 

To  provide  an  illustration  for  the  nature  of  the  operational  problem  and 
t.he  type  of  solution  sought,  consider  the  mission  trajectory  of 
Figure  21.  The  "A",  "B",  and  "C"  threat  categories  describe  three 

types  of  threats.  The  end  of  the  dotted  line  shows  the  current  loca¬ 
tion  of  the  aircraft.  At  this  point,  a  pop-up  threat  A2  is  detected. 
The  algorithm  generates  the  solid  trajectory  in  real  time  to  accomodate 
this  event,  Note  that  besides  generating  the  horizontal  (shown  in  the 
figure)  and  vertical  (not  shown)  plane  trajectories,  the  algorithm  also 
suggests  an  expendable  and  nonexpendable  resource  deployment  sche¬ 
dule,  depicted  in  the  lower  portion  of  t.he  figure.  On  the  upper  right 
the  fuel  used  is  shown;  the  probability  of  survival,  and  the  probability 
of  achieving  objectives  1  and  2  are  also  shown .  The  algorithms  devel¬ 
oped  by  SCI  (Vt)  email  the  use  of  dynamic  programming  methods  impel- 
mented  on  an  array  processor. 5  Modifications  to  these  algorithms  would 
be  required  if  a  serial  processor  is  used.  However,  Figure  21  does 
show  the  nature  of  the  output  that  ought  to  be  generated  by  a  meaning¬ 
ful  operational  algorithm. 

3.3.4  Threat  Modeling  Survey 

This  section  surveys  some  state-of-the-art  techniques  of  threat 
modeling.  Generally  speaking,  no  validated  threat  models  are  available 
at  this  time.  However,  it  <s  useful  to  review  the  modeling  methods  cur¬ 
rently  used.  As  noted  earlier  the  threat  must  be  described  in  terms 
of  spatial  elements  wherein  the  probability  of  surviving  transition  from 
one  cell  to  the  next  is  specified  as  a  function  of  aircraft  position  and 
velocity  with  respect  to  the  threat;  either  tabular  storage  or  curves 
could  be  used  for  this. 

One  approach  to  model  development  for  the  slow- allocation  mode  for  a 
SAM  site  is  to  perform  off-line  Monte  Carlo  runs  simulating  fly-past  at 
different  offsets  and  relative  heading  using  available  SAM  data  on  each 
Monte  Carlo  run.  The  A/C  cross  s.  ction  (radar/optical)  history,  F.CM 
effectiveness,  and  missile  effectiveness  must  be  chosen  independently. 
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In  this  way  the  unconditioned  probability  of  surviving  can  be  quantified 
probabilistically.  Models  generated  by  these  techniques  are  not  avail¬ 
able  at  this  time  and  are  the  subject  of  on-going  research  studies.  In 
lieu  of  these  validated  models,  a  stochastic  threat  model  can  be  devel¬ 
oped  as  detailed  in  Appendix  B .  It  is  suggested  that  this  stochastic 
model  might  be  a  satisfactory  means  of  tuning  up  the  threat-avoidance 
algorithm  by  using  a  range  of  numerical  values  for  the  threat  lethality 
parameters . 

3.3.5  Conclusions 


Threat  avoidance  and  threat  modeling  literature,  particularly  SCI 
work,  was  surveyed  to  determine  the  following: 

a.  No  validated  threat- survival  models  exists. 

b.  A  generalized  stochastic  threat  model  has  been  formulated 
(Appendix  B). 

c .  Operational  threat-avoidance  algorithms  should  provide  re¬ 
source  allocation  capability. 

d.  A  dynamic  programming  approach  on  an  array  processor  can 
be  used  to  implement  a  real-time  threat-avoidance  algorithm. 

Finally ,  it  is  noted  that  a  more  detailed  analysis  of  threat  avoidance  is 
beyond  the  scope  of  this  advisory  study  and  substantial  additional  work 
remains  to  be  performed  before  a  realistic  modification  of  the  existing 
trajectory  generator  algorithms  can  be  implemented. 

3.4  SIMULATOR  DEVELOPMENT 

The  major  areas  cf  simulator  development  were  the  following. 

a.  Applying  the  developed  six-degree-of-frsedom  KC-135  aero¬ 
dynamics  and  the  Stability  Augmentation  System  models  to  the 
closed  loop  .simulation . 

b.  Converting  the  ’ntegrated  Flight  Trajectory  Control  (1FTC) 
program  from  the  IBM-370  computer  to  a  PDP-11/70. 

c.  Modifying  the  control-display  hardware  for  the  transport 
cockpit. 

The  latter  two  tasks  were  accomplished  under  in-house  funding. 
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3.4.1  Aerodynamics  and  Stability  Augmentation  System 


The  F-4  aerodynamics  in  the  existing  IFTC  closed-loop  hybrid 
simulation  were  replaced  with  the  KC-135  aerodynamics. 

The  Stability  Augmentation  System  (SAS)  was  developed  and  checked 
out  with  the  aircraft  model.  The  mechanization  and  characteristics  of 
the  SAS  are  similar  to  existing  functions  of  the  digital  flight  control 
system  presently  installed  on  the  Speckled  Trout,6 

3.4. 1.1  SAS  Design 

The  SAS  structure,  as  defined  in  Boeing  Aircraft  Company  re¬ 
port,  "KC/E/RC-135  System  Characteristics  Missions  Simulator",  consists 
of  a  pitch  axis,  rcll  axis,  and  yaw  damper. 

The  pitch  axis  feedbacks  are  pitch  attitude  and  washed-out  pitch  rate. 
The  lateral  axis  feedbacks  are  roll  attitude  and  roll  rate.  The  yaw 
damper  feedback  is  washed-out  yaw  rate.  Gains  for  these  loops  are  ob¬ 
tained  using  a  standard  root  locus  analysis  program  that  is  resident  in 
the  PDP- 11/70  laboratory  computer. 

The  pitch  axis  block  diagram  is  shown  in  Figure  22. 


SAS  PITCH  AXIS  BLOCK  DIAGRAM 
FIGURE  22 


6Speckled  Trout  Software  Documents,  Specification  #5710,  Sperry 
Flight  Ssytems,  Phoenix,  Arizona. 
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In  this  diagram 

Ke  = 

pitch  attitude  feedback  gain 

K 

q 

pitch  rate  feedback  gain 

Q  = 

dynamic  pressure 

a_ 

6 

e 

(s)  = 

aircraft  short  period  mode  transfer 
function 

Selection  of  K„  and  K  depends  on  the  flight  condition.  The  block  -7c- 

0  V  ^ 

was  inserted  into  the  pitch  axis  control  in  order  to  adapt  the  gains  to 
varying  flight  conditions.  The  pitch  SAS  will  be  evaluated  for  both  a 
high  Q  and  low  Q  flight  condition.  The  flight  conditions  selected  are 
summarized  in  Table  IV. 


TABLE  IV 

SUMMARY  OF  A/C  FLIGHT  CONDITIONS 


Aircraft 

Flight 

Condition 

Parameter 

Low  Q 

High  Q 

h 

5000' 

25000' 

M 

.25 

.75 

P 

.002049 

.001065 

V 

274.8  ft/sec 

763.4  ft/sec 

Q 

77.4  lbs/ft2 

310.3  Ib/sec 

a . 

trim 

7.19° 

-.36° 

SWCP, 

trim 

-7.13° 

-2.27° 

T.  • 

trim 

7711  lbs 

12728  lbs 
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The  short  period  mode  transfer  functions  are  derived  in  Appendix  A. 
For  the  low  Q  flight  condition: 


6~  (s) 
e 


-  . 921 (s  +  .717) 

P  +  1,75s  +2.05 


(99) 


For  the  high  Q  flight  condition : 


SL.  f'o'i 
6  '  j 
e 


- 2 . 6 (s  +  .83) 

Z*  +  2.45  +  8.17 


(100) 


The  root  loci  for  both  flight  conditions  were  obtained  by  breaking  the 
innermost  loop  and  are  shown  in  Figure  23  for 


K 


0 
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Selecting  =  1  results  in  good  transient  response  characteristics, 
The  characteristics  of  the  dominant  complex  poles  are  summarized  in 
Table  V. 


TABLE  V 

FITCH  SAS  RESPONSE  CHARACTERISTICS 


Flight  Condition  Damping  Ratio  (£) 


Low  Q 
High  Q 


.79 

.49 


Natural  Frequency  (uj^) 

3.15  rad/sec 
3.47  rad/ sec 
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The  yaw  damper  block  diagram  is  shown  in  Figure  24. 


SAS  YAW  DAMPER  BLOCK  DIAGRAM 
FIGURE  24 


Derivation  of  the  aircraft  transfer  function  (s)  is  shown  in 
Appendix  A  °R 

For  the  low  Q  flight  condition 

r_  ,  v  _  - . 447 (s+1 . 25 ) [ (s- ■ 0934) 3  +  (.437)2) 

6r  "  (s+1. 26)  (s+. 00589)  [(b+.00586P  +  (. 954)"^ 

For  the  high  Q  flight  condition 

r_  ,  .  _  -1.54(s+1.44)[(s-.0976)2+(,  436) 2  ] 

6r  (S}  ~  ( B+ 1 . 49 ) ( s+ . 00668 ( ( ( s+709 1 3 ) 2+ ( 1 . 69 ) *7 

The  root  loci  for  both  conditions  are  shown  in  Figure  25.  KR  =  2  was 
chosen  for  maximum  damping  of  the  dominant  complex  poles. 


The  roll  axis  SAS  is  shown  in  block  diagram  form  in  Figure  26. 


les 


Servo  A/C 


ROLL  AXIS  BLOCK  DIAGRAM 
FIGURE  26 


In  Figure  26, 


K,  =  roll  attitude  feedback  gain 
<)> 


K 


roll  rate  feedback  gain 


*-<»> 

a 


r-+6. 


is  the  roll  attitude  to  aileron  transfer  function  with  the  yaw  loop 
closed  and  is  derived  in  Appendix  A.  The  root  locus  for  the  lateral 
SAS  was  obtained  by  breaking  the  innermost  loop  and  setting 


K. 

J& 

K 

P 


1.5 


The  root  loci  for  both  flight  conditions  are  shown  in  Figure  27.  A 

value  of  k  =  2.5  was  chosen  for  the  roll  rate  gain.  The  characteris- 
P 

tics  of  the  cominant  complex  poles  are  summarized  in  Table  VI. 


TABLE  VI 

LATERAL  SAS  RESPONSE  CHARACTERISTICS 


Flight  Condition 

Damping  Ratio  (£) 

Natural  Frequency  (u)n) 

Low  Q 

.4 

.695  rad/sec 

High  Q 

.2 

1.53  rad/sec 

33 
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3.4. 1.2  SAS  Mechanization 


The  SAS  dynamics  were  simulated  on  the  laboratory  analog 
computer  with  several  of  the  gains  computed  digitally  and  set  on  the 
appropriate  pots.  The  SAS  block  diagrams  written  in  equation  form  be¬ 
come 


q  ^s+l^s+4.55^  f  ^eCMD+6^  K. 


q  J 


Q  q 


'_L_i t  15.7  s  /100-. 
‘s+lns+15.7n  Q  J 


K, 


^CMD-^  K  "  p 


(IT*  kp  ( 


10 ) 
s+10"' 

= 

'6 

e 

(103) 

10  ) 
svl0j 

- 

6r 

(104) 

10 ) 
s+10; 

= 

6 

a 

005) 

Definitions  of  the  terms  used  in  equations  99  through  105  are  given  in 
Table  VII. 


These  equations  are  scaled  by  redefining  the  problem  variables  as  the 
ratios  of  the  original  variables  divided  by  their  maximum  magnitude. 
For  a  100-volt  analog  machine,  a  scaled  variable  with  value  1  is  equi¬ 
valent  to  100  volts.  Maximum  magnitudes  assigned  to  the  problem  vari¬ 
ables  are  shown  in  Table  VIII. 


Rewriting  the  SAS  equations 

A  «;>;  100K  in 

l1-835  ciTT>ci+trl5)  qt  +  2-3375  r  (We)!  tj 

q 


ov  ri  /  s  s /  15.7  ■.  /100-, r  10  i 
,8KR  (s+r  (s+15.7)(  Q  )ts+10) 


(106) 


'V  W*1  •  V  (T)  2'865  ‘iiTo1  = 

The  analog  implementation  of  these  equations  is  shown  in  Figure  28. 


3.4. 1.3  SAS  Performance 


The  SAS  design  and  mechanization  was  verified  by  recording 
the  longitudinal  and  lateral  SAS  responses  to  pulse  elevator  and  aileron 
deflections,  respectively.  Figures  29  and  30  show  these  responses, 
Table  IX  compares  the  recordings  with  the  calculated  SAS  transient  re¬ 
sponse  characteristics . 
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TABLE  VII 

TERM  DEFINITION,  EQUATIONS  99  through  105 


pitch  attitude  feedback  gain  (rad/rad) 
pitch  rate  feedback  gain  (rad/rad/sec) 
yaw  rate  feedback  gain  (rad/rad/sec) 
roll  attitude  feedback  gain  (rad/rad/sec) 
roll  rate  feedback  gain  (rad/rad/sec) 
elevator  deflection  angle  (rad) 
rudder  deflection  angle  (rad) 
aileron  deflection  angle  (rad) 
laplace  transform  operator 
dynamic  pressure  (lbs/ft2) 
altitude  (ft) 

Mach  number 

air  density  (slugs/ft3) 

total  linear  velocity  (ft/sec) 

angle  of  attack  (deg) 

stabilizer  deflection  angle  (deg) 

engine  thrust  (lbs) 

yaw  rate  (rad/sec) 

roll  rate  (rad/sec) 

pitch  rate  (rao/sec) 

pitch  attitude  (rad) 

damping  ratio 


TABLE  VII 

TERM  DEFINITION,  EQUATIONS  99  through  105  (Continued) 


U) 

n 

- 

natural  frequency  (rad/sec) 

c 

= 

wing  chord  (20.16  ft/) 

S 

= 

wing  reference  area  (24.33  ft2) 

IYY 

- 

pitch  moment  of  inertia  (slug-ft2) 

m 

- 

aircraft  mass  (4204  slugs) 

g 

- 

gravity  constant  (32.2  ft/sec2) 

b 

- 

wing  span  (130.83  ft) 

[xx 

- 

roll  moment  of  inertial  (1.985xl06  slug-ft2) 

*ZZ 

= 

yaw  momemt  of  inertia  (4.608x10®  slug-ft2) 

:xz 

n 

product  of  inertia  (.147x10®  slug-ft2) 

TABLE  VIII.  PROBLEM  VARIABLE  SCALING 
Symbol  Maximum  Value  Scaled  Variable 


.785  rad/sec 

q’ 

= 

[.785j 

.4278  rad 

6  ' 
e 

'  6e 

.4278 

.349  rad/sec 

r' 

= 

L . 349 J 

436 

6R 

= 

~  «r1 

.  436_ 

"  “ 

.349  rad 

6  ' 
a 

= 

.349 
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STABILITY  AUGMENTATION  SYSTEM  CIRCUIT  DIAGRAMS 

FIGURE  28 


TABLE  IX 

SAS  DESIGN  VERIFICATION 


Axis 

Flight 

Condition 

Measured  SAS 
Response 

Calculated  SAS 
Response 

4  u>  (rad/sec) 

4 

u)n  (rad/sec) 

Longitudinal 

Low  Q 

.8  3.14 

.79 

3.15 

High  Q 

.5  3.3 

.49 

3.47 

Lateral/Directional 

Low  Q 

.4  .698 

.4 

.695 

High  Q 

.14  1.57 

.2 

1.53 

3. 4. 1.4  KC-135  Aerodynamics 

Six-degree-of-freedom  aircraft  equations  of  motion  were  de¬ 
rived  to  simulate  the  KC-135  aircraft  aerodynamics. 


The  equations  of  motion  of  an  airframe  referred  to  Eulerian  axes  are: 


• 

U 

1 

tj 

X 

_ i 

- 1 

O 

n 

i 

•o 

m  m 

u 

• 

V 

1 

ZFV 

+ 

-r  0  p 

V 

mass 

Y 

w 

-  m 

1 

ft 

_ 1 

q  -p  0 

»  m 

w 

m  • 
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P 

1vl° 

IL 

0 

r 

-q 

1XX° 

-Ixz 

p 

♦ 

q 

= 

0  XYY 

0 

IM 

+ 

-r 

0 

p 

0  ^ 

0 

q 

• 

r 

•  ■ 

."^XZ0 

HH 

N 

N 

1 _ 

ZN 

q 

“P 

0 

m 

.^xz0 

llZ  . 

r 
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The  Euler  angle  rates 


$ 

"1  sin$tan0  cos<f>tanG 

’p" 

e 

= 

0  cos$  -sin$ 

q 

_0  sin$/cos0  cos$/cos0_ 

_  r  _ 

(108) 


(109) 
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The  external  forces  and  moments  acting  on  the  aircraft  consist  of  aero¬ 
dynamic,  propulsion,  and  gravity  components.  The  translational  equa¬ 
tions  are  solved  using  the  wind  axis  system  since  it  is  well  known  that 
much  lower  accuracy  and  speed  demands  are  made  on  computer  mechani¬ 
zation  than  a  body  axis  solution,7  The  aerodynamic  forces  defined  in 
the  stability  axis  system  and  che  gravity  and  propulsion  forces  defined 
in  the  body  axis  system  are  transformed  into  the  wind  axis  coordinate 
system  as  shown  below. 


X 

w 

^x 

cosp 

sinp  0 

X 

s 

cosa 

0 

sina 

T 

-sin6 

Y 

w 

=: 

IFy 

- 

-sinp 

cosP  0 

Y 

s 

+ 

0 

1 

0 

0 

+w 

cos0sin<t> 

Z 

w 

0 

0  1 

Z 

s 

-sinof 

0 

cosa 

0 

COS0COS<J> 

By  definition,  the  velocity  vector  of  the  aircraft  (V^,)  is  along  the 
axis  and,  therefore,  Ww=Ww=Vw=Vw=0.  We  have  then 


mVn 


mV„  r 
T  w 


mVT  qw 


or  r  =  Y  /mVT 
w  w  w  T 


(111) 


or  q  =  -Z  /mV-, 
w  v  w  T 


where  rw  and  qw  are  flight  path  axis  angular  rates  about  the  Zw  and 
Xw  axes,  respectively.  Equations  for  c  and  p  are  derived7  from  the 
geometrical  relationships  between  axis  systems  (Figure  31). 


feasibility  Study  for  Integrated  Flight  Trajectory  Control  (Airlift). 
LSI  Report  #  GRR-008-0177A,  AFFDL-TR-77- 120 ,  contract  F33615-74-C- 
3083 . 
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body  stability,  and  wind  axis  system  geometries 

’  FIGURE  31 


«  *  |<v08fi '  V^i  t  s;| 

Yw 

P  ~  mV  ’  rs 
s 


(112) 


After  expanding  and  simplying,  the  following  translational  equations 
result: 


P 

a 


V" 

Vvt  - 

Zw/(mVTcosP)  +  qg  -  Pfitanp 


(113) 


! 


The  rotational  equations  of  motion  are  computed  in  the  stability  axis. 
After  simplification  we  have 


(Ilrb  +  l^%  +  + 

I5rbPb  +  Ie(rJ-pj)  +  I7Mb 

(IsPb  +  Iqrb)qb  +  l4Lb  +  IloNb 


(114) 


where 
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The  aerodynamic  forces  and  moments  are  defined  below. 


-^vt2scd 


SspVTzSbC£ 


JspvT2scY 


(116) 


-JsPvT2scL 


JspVTzSbCn 


The  dimensionless  aerodynamic  coefficients  (Cn,  Cv,  CT  .  C0,  C  ,  C  ,) 

u  i  l  x  in  n 

were  obtained  with  assumptions  defined  in  Section  3.2.1. 

With  reference  to  Appendix  A,  the  Aerodynamic  coefficient  equations  for 
simulation  requirements,  reduce  to  the  following; 


1  -  .—i 
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The  terms  of  the  aerodynamic  coefficient  eolations  are  stored  in  compu¬ 
ter  memory  in  a  table  lookup  format.  Term  definitions  have  been  given 
previously  after  equation  59. 


'1  KC-  135  aerodynamics  were  incorporated  in  a  complete  man-in-the- 
loop  simulation  using  a  hybrid  analog/digital  computer  and  our  two-man 
cockpit. 

The  driving  function  inputs  to  the  simulation,  as  mechanized  in  the 
block  diagram  of  figure  32,  are  the  cockpit  stick  and  throttle  controls 
manipulated  by  the  "pilot" . 

The  computations  performed  by  the  blocks  indicated  with  an  asterick 
have  been  previously  described.  The  remaining  blocks  perform  the  fol- 
U  ng  \>iaputations : 

n  Body  Axis  Velocity  Computations 


U  =  •  eosacosjl 

V  =  VT  •  sinp 

W  =  VT  ♦  sinacosp 

□  Euler  angles  and  direction  cosine  matrix  computations 
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□  Inertial  velocity  computations 
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Aircraft  position  computations 
x  - 


y 

h 


-rVEgdt 

J-VDdc 


(122) 


(123) 


(124) 


!  2 


98 


ij 

If 


□  Air  data  computations 
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Control  surface  computations 
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The  translantional  and  rotational  equations  of  motion  and  the  control 
surface  dynamics  programmed  on  the  analog  computer  were  scaled  by 
redefining  the  problem  variables  as  the  ratios  of  the  original  variables 
divided  by  their  maximum  magnitude.  For  a  100-volt  analog  machine,  a 
scaled  variable  with  value  1  is  equivalent  to  100  volts.  Maximum  magni¬ 
tudes  assigned  to  the  problem  variables  are  shown  in  Table  X. 

Solving  the  rotational  equations  of  motion  using  the  redefined  problem 
variables  we  have 

p*  =  .2123  r'q'  +  .0320  p  'q '  +  .6433  Lg  +  .0205  Ng 

q'  =  .4887  r'p'  +  -0122  t'2  -  .0616  p  '2  +  .6802  Mg  (128) 

r'  =  -.2075  p'q'  “  •  0320  r'q'  +  .0461  Lg  +  .6232  N' 

The  translational  equations  of  motion  become 
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and  the  control  surface  dynamics  are 


6 

c 

- 

(6' 

eCMD 

-6^/.l 

-1.  5 

6' 

e 

§  .6122 

6' 

a 

= 

(6a 

cMr 

-6a')/.! 

-1.  s 

5' 

a 

S  1. 

(130) 

= 

(6' 

CMD 

~6r)/a 

-1.  § 

6r 

«  1. 

-.0357 

5  SWCP  s  1 

JW 


VTcosP 


+  q '  -  1 .4981  p '  tanP 


100 


. ^  a  -1r 


H=t*r 


SYMBOL 


MAXIMUM  VALUE 
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The  preceding  equations,  together  with  the  cockpit  trunking,  were  im¬ 
plemented  as  shown  in  the  analog  circuit  diagrams  contained  in 
Appendix  C. 

The  digital  computations  are  performed  every  50  milliseconds.  A  func¬ 
tional  flowchart  of  these  computations  is  contained  in  Appendix  C . 

Aircraft  response  to  pulse  elevator  and  pulse  aileron  deflections  are 
shown  in  Figures  33  and  34  for  both  a  high  Q  and  a  low  Q  flight  condi¬ 
tion.  Table  XI  compares  the  simulated  and  desired  aircraft  models  by 
presenting  the  characteristics  of  the  dominant  lateral  and  longitudinal 
aircraft  modes. 


TABLE  XI 

AIRCRAFT  DYNAMICS  VERIFICATION 


3.4.2  Rehosting  of  IFTC  Algorithms 


The  IFTC  algorithms,  developed  on  previous  AFFDL  programs, 
were  incorporated  in  a  closed-loop  simulation  through  a  general  purpose 
IBM-370  computer.  This  IBM  computer  was  not  dedicated  to  the  Simula¬ 
tion  Laboratory  but  was  accessed  through  a  real-time  monitor. 


The  KC-135  cockpit  simulator  was  upgraded  by  rehosting  the  IFTC  algo¬ 
rithms  from  the  IBM-370  computer  to  a  dedicated  PDP-11/70  simulation 
computer. 
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Because  of  the  unreliability  of  the  real-time  monitor,  this  change  in  the 
mechanization  of  the  simulation  increased  the  time  availability  and  effi¬ 
ciency  of  the  simulator.  Simultaneously  with  this  rehosting  of  the  soft¬ 
ware,  the  algorithms  were  upgraded  to  include  racetrack  pattern  and 
spiral  climb  capability  of  the  aircraft. 

The  planned  incorporation  of  the  threat  avoidance  algorithms  was  never 
completed  because  the  specific  form  of  the  algorithms,  suitable  for  im¬ 
plementation  on  our  airborne  avionics  computer,  was  not  available. 

The  review  of  the  software  during  the  rehosting  process  provided  sev¬ 
eral  areas  for  simplification  and  improvement  of  the  software  to  improve 
portability  and  performance. 

3.4.3  Controller/Display  Modification 

The  controller/display  in  the  cockpit  was  modified  by  addition  of 
a  refresh  memory  to  provide  increased  performance  of  the  display  and 
to  provide  modified  row/column  matrix  select  keys  to  improve  the  flexi¬ 
bility  of  menu  select  and  modification  functions. 

Figure  35  is  a  faceplate  of  the  controller/display  currently  installed  in 
the  transport  cockpit. 
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RECOMMENDED  FLIGHT  TEST  PLAN 


Although  the  optimal  flight  trajectory  control  algorithms  are  not  yet 
state-of-the-art  for  systems,  the  baseline  classical  integrated  flight  tra¬ 
jectory  control  algorithms  have  been  simulated  but  not  demonstrated  in 
an  operational  aircraft.  The  purpose  of  this  flight  test  plan  is  to  de¬ 
fine  a  program  that  will  use  the  Speckled  Trout  aircraft  to  demonstrate 
the  operational  utility  of  the  conventional  IFTC  algorithms  in  the  KC-135 
tanker  mission. 

4.1  FLIGHT  TEST  PLAN  OVERVIEW 

The  major  tasks  and  other  relationships  for  the  flight  test  plan  are 
shown  in  Figure  36,  Task  Flow  Diagram. 

Typical  KC-135  Tanker  mission  scenarios  are  analyzed  for  operational 
requirements  of  the  avionics  system.  These  requirements  are  then  com¬ 
pared  to  the  unique  capabilities  of  the  IFTC  algorithms  to  establish  the 
principal  flight  test  requirements  of  time  control  and  in-flight  mission 
redirect.  The  avionics  system  mechanization  and  flight  test  scenarios 
that  are  necessary  to  demonstrate  the  test  requirements  are  limited  to 
mechanical  and  operational  constraints  of  the  flight  test  aircraft.  The 
mechanized  constraints  are  primarily  cockpit  panel  space  and  rack  avail¬ 
ability  for  the  equipment.  The  operational  constraints  are  the  limited 
time  availability  for  modifications  to  the  aircraft  and  the  requirement  to 
maintain  a  flightworthy  operational  avionics  suite  throughout  this  test 
program.  Our  basic  concept  to  achieve  these  difficult  operational  re¬ 
strictions  is  to  configure  a  three-phase  modification  and  test  procedure. 

Air  Force  personnel  will  perform  the  actual  equipment  installation,  check 
out  the  data  collection,  perform  the  data  and  analysis  with  support  from 
LSI .  The  Air  Force  is  currently  developing  the  software  for  an  opera¬ 
tional  flight  program  and  critical  elements  of  the  avionics  update  pro¬ 
gram.  This  system--which  will  be  tested  functionally  on  the  bench  and 
operationally  in  a  fixed-base  simulator--will  contain  the  IFTC  algorithm 
functional  capabilities.  This  system,  or  an  existing  LSI  commercial- 
hardware  equivalent  with  software  modifications,  is  proposed  for  this 
flight  test. 
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4.2  FLIGHT  TEST  FUNCTIONAL  REQUIREMENTS  ANALYSIS 


The  mission  of  the  KC-135  aircraft  is  to  provide  air  refueling  and/ 
or  air/land  stores  delivery  anywhere  in  the  wond  while  operating  with  a 
three-man  crew.  This  general  mission  requirement,  and  the  specific 
operational  scnenarios  of  three  typical  missions8,  Doth  explicitly  and  im¬ 
plicitly  define  the  functional  requirements  for  the  avionics  system  of  the 
aircraft.  The  detailed  functional  requirements  for  navigation,  guidance, 
flight  plan  management,  communications,  performance  management,  con¬ 
trol  and  display  are  outlined  in  the  following  paragraphs.  These  de¬ 
tailed  functional  requirements  provide  the  basis  for  the  flight  test  sys¬ 
tem  mechanization  (Section  4.3). 


4.2.1  General  Navigation 


A  fundamental  avionics  system  requirement  is  the  capability  to 
determine  the  aircraft  position  and  velocity,  with  respect  to  the  current 
mission  plan  and  to  other  aircrrft,  with  a  minimum  crew  workload. 


These  general  capabilities  imply  the  following  specific  requirements: 

a.  Communications  with  mission  control  and  ower  mission  air¬ 
craft. 


fc  Measurement  of  relative  position/velocity  vectors  with  respect 
to  other  mission  aircraft  using  cooperative  and/or  non-cooper¬ 
ative  sensors. 


c.  Knowledge  of,  and  ready  access  to,  current  navigation  data 
associated  with  foreign  and  domestic  route  structures ,  air¬ 
ways,  waypoints,  navigation  aids  and  terminal  areas  for  mis¬ 
sion  planning, 

d.  Utilization  of  existing  and  future  civil  and  military  external 
navigation  aids. 


Because  enemy  action  may  shut  down  these  external  navigation  aids  dur¬ 
ing  all  or  part  of  some  missions  (e.g.,  Mildenhall  EWO  mission)  the  air¬ 
craft  must  achieve  the  minimum  acceptable  navigation  performance  with 
only  self-contained  dead  reckoning  equipment. 

For  minimum  crew  workload  the  navigation  system  must  have  the  capa¬ 
bility  to  automatically  select  the  best  available  combination  of  internal/ 
external  navigation  sensors,  throughout  the  mission,  and  integrate  their 
outputs  to  obtain  a  best  estimate  of  the  aircraft  navigation  state. 


8R.P.  Madero,  TAACE  Scenario  Update,  Bunker-Ramo  Document 
No.  4506-020-5100-9,  May  1979. 


4.2.2  Guidance 


The  basic  air  refueling  mission,  particularly  those  involving  SAC 
receivers,  requires  precise  ARCP/ARCT  rendezvous.  Explicit  in  this 
basic  mission  requirement  ere  functional  avionics  system  capabilities  for 
precise  mission  time  control  and  in-flight  mission  redirect. 

The  general  requirement  for  precise  time  control  implies  the  following 
specific  capabilities  for  the  avionics  system: 

a.  Automatic  generation  of  a  reference  flight  path  in  four  dimen¬ 
sions  that  is  constrained  only  by  the  dynamics  of  the  aircraft 
and  the  waypoints  specified  by  missions  planning. 

b.  Automatic  prediction  of  the  future  position  and  velocity  of  the 
aircraft  from  the  estimates  of  the  present  aircraft  state  as 
generated  by  the  navigation  system. 

c.  Automatic  generation  of  steering  and  throttle  commands  that 
will  minimize  errors  between  the  aircraft  present  position  and 
the  desired  position  of  the  mission  plan. 

The  general  requirements  for  precise  in-flight  mission  redirect  implies 
specific  capabilities  for  easy  reference  path  modification  by  manual  in¬ 
puts  from  the  crew  (see  Controller/Display,  Section  4.2.5).  These 
changes  in  reference  path  may  be  required  to  avoid  unforeseen  threats 
such  as  weather  (mildenhall  EWO  mission)  with  minimum  disruption  of 
the  SAC  receiver  mission  plan.  Or,  the  reference  path  may  require  a 
change  in  order  to  reach  out  in  performing  a  rendezvous  in  minimum 
time  with  fighter  aircraft  returning  from  the  battle  area  with  unsuffi¬ 
cient  fuel  to  make  the  original  ARCP  (BODO  contingency  mission). 

Also  implicit  in  the  general  requirements  for  precise  time  control  and 
in-flight  mission  redirect  is  the  specific  requirement  fo  minimize  tanker 
fuel  consumption  in  order  to  maximize  fuel  available  for  offloading  (see 
Performance,  Section  4.2.4).  This  specific  requirement  results  from  the 
fact  that  the  KC-135,  in  flight,  draws  fuel  from  the  same  fuel  supply 
that  is  used  to  refuel  the  receivers. 


4.2.3  Flight  Plan  Management 


The  functional  requirements  of  the  navigation  system  for  know¬ 
ledge  of  and  access  to  a  current  worldwide  navigation  data  base  of  the 
guidance  system  to  generate  a  reference  flight  path  through  mission 
waypoints  specified  in  four  dimensions  imply  the  following  two  functions 
for  flight  plan  management: 


a. 


Accept,  edit,  and  store  the  navigation  data  associated  with 
commercial  and  military  route  structures,  airways,  waypoints, 
navigation  aids  and  terminal  areas.  These  data  must  be  up¬ 
dated  every  month. 

b.  Accept,  edit  and  store  the  pai  dcular  mission  waypoint  se¬ 
quence  at  the  initiation  of  each  mission. 

For  minimum  crew  workload  and  aircraft  system  reaction  time  (preflight 
preparation),  this  data  should  be  loaded  automatically  into  the  avionics 
system  (see  Controller/Display) . 

4.2.4  Performance 

The  KC-135  aircraft,  in  flight,  draws  from  the  same  fuel  supply 
that  is  used  to  refuel  the  receivers.  The  functional  requirement  of  air¬ 
craft  performance  management  is  to  minimize  the  fuel  consumed  by  the 
tanker  in  order  to  maximize  the  fuel  available  for  offloading.  In  order 
to  meet  this  functional  requirement,  t,  »  mission  time  constraint  of  ARCT 
should  be  compatible  with  the  mission  position  constraint  of  ARCP.  For 
example,  at  250  KIAS  refueling  velocity,  a  time  constraint  of  ±  1  minute 
corresponds  to  a  position  constraint  on  ARCP  of  ±  5  nautical  miles. 
The  associated  fuel  savings  is  approximately  200  pounds  per  minute. 

Minimum  fuel  consumption  enroute  implies  the  following  required  capa¬ 
bilities  : 

a.  Control  of  the  fuel  offloaded  from  various  tanks  in  order  to 
maintain  proper  c.g.  of  the  aircraft  to  minimize  aerodynamic 
drag  of  the  aircraft. 

b.  Flight  mode  control,  within  constraints  of  the  mission  planning 
to  achieve  a  specified  aRCP/ARCT,  for  minimum  fuel  con¬ 
sumption.  These  modes,  for  example,  include  maximum 
economy  climb-out,  cruise,  and  descent,  and  wind  trades  to 
obtain  operate  ig  altitude  for  minimum  fuel  burn  during  cruise, 

c.  Operation  of  the  engines  at  optimum  engine  pressure  ratio  for 
minimum  fuel  consumption  through  all  flight  regimes. 

In  addition,  in  order  to  assure  maximum  time  in  the  rendezvous  refuel¬ 
ing  pattern  (BODO  contingency  mission),  the  Aircraft  Performance  Man¬ 
agement  System  must  continously  compute  projected  fuel  remaining  at 
end  of  mission. 


4.2.5  Controller/Display 


The  principal  function  of  Controller/Display  is  to  allow  efficient 
management  of  the  aircraft  by  the  crew.  In  general,  this  requires  both 
specific  dedicated  displays  and  integrated  controller/displays  that  are 
functionally  interactive  with  the  crew.  Specifically,  these  requirements 
include : 

a.  Integrated  alphanumeric  controller/ displays  for  monitoring  and 
controlling  system  operational  modes,  flight  plans,  and  air¬ 
craft  system  status. 

b.  Integrated  graphics  displays  for  showing  the  aircraft  situation 
in  both  horizontal  and  vertical  planes  and  in  relation  to  the 
current  flight  plan,  identified  threats  such  as  weather,  and 
friendly  and  bogey  aircraft  in  the  operational  area. 

c.  Dedicated  displays  for  warning  of  critical  system  malfunction 
and  for  monitoring  and  control  of  aircraft  dynamics  and 
engine  performance  (attitude,  altitude,  airspeed,  engine  pres¬ 
sure  ratio). 

d.  Dedicated  controls  for  automatic  entry  of  initial  flight  plan 
data  and  navigation  data  base. 

4.2.6  Communications 


The  aircraft  must  have  direct  communications  with  other  opera¬ 
tional  units  in  the  mission,  the  fuel  receivers,  and  ground  control. 
The  required  information  includes  status  of  area  friendly  and  enemy  air¬ 
craft,  weather  patterns  that  may  affect  the  mission,  and  ground  control 
missions  redirect  information. 

4.3  FLIGHT  TEST  SYSTEM  MECHANIZATION 

The  proposed  flight  test  avionics  system  mechanization  is  based  on 
our  commercial  Performance  Navigation  Computer  System  (PNCS),  devel¬ 
oped  for  the  Boeing  Aircraft  Company  Models  727  and  737  aircraft  and/ 
or  the  system  currently  being  developed  by  ASP/EN  for  the  KC-135  avi¬ 
onics  update  program. 

4.3.1  Functional  Description 

The  PNCS  shown  in  Figure  37  is  a  functionally  integrated  per¬ 
formance/guidance  mission  management  system.  The  system  software 
that  is  currently  being  adapted  and  modified  for  the  KC-135  avionics 
update  program  includes  the  trajectory  generator,  path  prediction,  and 
guidance  algorithms  of  the  Integrated  Flight  Trajectory  Control  system. 
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The  functional  block  diagram  of  the  proposed  flight  test  system  is 
shown  in  Figure  38. 


FUNCTIONAL  BLOCK  DIAGRAM,  FLIGHT  TEST  SYSTEM 

FIGURE  38 

The  cockpit  management  function  allows  pilot  interactive  interface  with 
the  mission  management  system  through  a  conventional  controller/display 
unit  and  a  multifunction  display  with  graphics  capability.  This  inter¬ 
face  allows  control  of  the  system  by  access,  insertion,  and  modification 
of  data  in  the  flight  plan  management,  guidance,  and  performance  func¬ 
tions  of  the  system. 

The  flight  plan  management  function  includes  flight  plan  selection,  edit¬ 
ing,  insertion  or  deletion,  and  waypoint  sequencing. 
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The  flight  plan  management  function  maintains  two  flight  plans:  en¬ 
gaged  and  temporary.  The  engaged  flight  plan  controls  the  aircraft 
through  guidance  and  cockpit  management  functions.  The  editing  func¬ 
tion  of  the  engaged  flight  plan  is  limited  to  providing  supplementary  in¬ 
formation  from  the  navigation  data  base  to  the  cockpit  management  func¬ 
tion  to  define  the  airports,  navaids,  and  waypoints  associated  with  the 
plan.  The  temporary  flight  plan  may  be  modified  by  waypoint  insertion 
and/or  deletion  through  interaction  with  the  cockpit  management  func¬ 
tion  and  using  data  from  the  navigation  data  base  and  guidance  func¬ 
tion.  The  temporary  flight  plan  may  become  the  engaged  flight  plan  or 
request  from  the  cockpit  management  function. 

The  waypoints  of  both  the  engaged  and  temporary  flight  plans  are  se¬ 
quenced,  as  a  function  of  the  current  state  of  the  aircraft,  on  command 
from  the  guidance  function.  The  major  guidance  functions  are  path 
prediction,  path  generation,  and  steering. 

Path  prediction  processing  generates  time-of-flight  parameters  based  on 
the  specified  waypoint  sequence  of  the  flight  plan  and  the  desired  per¬ 
formance  mode  of  operation. 

Guidance  path  processing  generates  the  reference  aircraft  flight  path  to 
correspond  to  the  waypoint  sequence  that  defines  the  flight  path.  This 
function  also  uses  the  current  aircraft  state  from  the  navigation  func¬ 
tion  to  determine  the  progress  of  the  aircraft  in  relation  to  the  desired 
or  reference  flight  path.  Guidance  path  processing  also  generates  a 
modification  to  the  refueling  holding  point  pattern  to  make  the  KC-135 
turn  out  ahead  of  the  receiver  aircraft  for  refueling  rendezvous,  and 
for  capture  of  the  flight  path  on  takeoff  and  ILS  beam  on  landing.  The 
steering  function  generates  lateral  and  longitudinal  steering  commands 
and  throttle  commands. 

The  navigation  function  integrates  aircraft  state  measurement  data  from 
on-board  dead- reckoning  and  ground-based  position  sensors.  The  navi¬ 
gation  function  draws  on  a  supply  of  reference  data  that  is  carried  on 
board.  This  navigation  data  base  includes  location  and  characteristics 
of  ground-based  navigation  aids,  standard  approach  and  departure 
routes  in  landing  sites,  and  standard  airways. 

The  performance  functions  are  takeoff,  c.g.  control,  flight  mode  con¬ 
trol,  and  performance  data  computation. 

Takeoff  calculates  data  to  predict  successful  takeoff  based  on  aircraft 
weight,  outside  air  temperature,  wind  vector,  runway  length  and  condi¬ 
tion  and  possible  loss  of  one  engine. 

The  c.g.  calculation  computes  and  controls  total  c.g.  of  the  aircraft 
based  on  fuel  burned  and  offloaded  from  each  fuel  tank. 


The  flight  mode  computes  target  EPR  and  airspeed  commands  for  all 
flight  modes:  Takeoff,  climb,  cruise,  descent,  hold,  maximum  continu¬ 
ous,  go-around,  and  turbulence. 

The  performance  processing  generates  data  used  by  the  flight  crew  to 
investigate  various  flight  profiles.  This  data,  which  would  normally  be 
obtained  by  reference  to  a  flight  operations  manual,  consists  of  data 
load,  flight  level  intercept,  flight  level  investigation,  ground  speed, 
range  and  endurance,  fuel,  temperature,  landing  flap/speed,  trip  plan¬ 
ning,  and  wind.  This  data  is  shown  to  the  pilot  either  in  the  form  of 
a-N  advisories  or  in  its  direct  effect  on  the  generated  flight  path  on 
the  map  display. 

4.3.2  Physical  Description 

The  phase  modular  character  of  the  proposed  flight  test  allows 
initial  limited  objective  tests  with  a  system  mechanization  that  provides 
less  than  the  full  functional  capability  required  for  the  total  mission 
management  system.  This  system  partitioning,  however,  also  provides  a 
piecemeal  approach  to  implementing  a  new  avionics  system  in  the  test 
aircraft  so  that  down-time  for  installation  and  checkout  are  limited  and 
flight  test  risks  are  reduced. 

Three  successively  more  complete  system  configurations  are  proposed 
for  flight  test:  baseline,  interim,  and  full-up  mission  management 
system . 

4.3.2. 1  Phase  One  Configuration:  Baseline  System 

For  Phase  One  of  the  flight  demonstration  program ,  a  single¬ 
thread,  baseline  system  will  be  installed,  This  baseline  system  consists 
of  the  following  units: 

□  One  Mission  Management  Computer  Unit  (MMC) 

□  One  Fuel  Savings  Advisory  System  Computer  (FSAS) 

□  One  Controller/Display  Unit  (CDU) 

The  inputs  required  from  the  auxiliary  equipment  currently  on  the 
Speckled  Trout  aircraft  to  generate  the  fuel  savings  (performance)  func¬ 
tions  are  the  following: 

a.  An  accurate  source  of  total  air  temperature  data  which  is 
available  from  the  Honeywell  CADC  (HG-180/W-747) . 
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b.  Air  data  inputs  of  pressure  altitude  from  the  CADC. 

c.  Calibrated  or  indicated  airspeed  from  the  CADC. 

d.  Engine  bleed  air  status  for  anti-ice  and  cabin  pressure. 

To  generate  the  guidance  and  navigation  functions  ,  the  required  inputs 
are  the  following: 

a.  Magnetic  heading  for  heading/airspeed  dead  reckoning  naviga¬ 
tion  . 

b.  Range  and  bearing  from  ground-based  reference  points  pro¬ 
vided  by  VOR/TACAN  stations, 

c.  Present  position  and  velocity  information  from  the  inertial  nav¬ 
igation  system  (INS). 

The  baseline  system  output  information  will  be  supplied  as  advisory 
flight  instrument  indications  to  the  copilot  only.  These  advisory  indica¬ 
tions  will  be  as  follows: 

a.  Copilot's  ADI  fast/slow  needle  or  MACH/1  AS-driven  bug  for 
time  control. 

b.  Copilot's  HS1  for  V-NAV  and  L-NAV  guidance  as  follows: 

□  The  course  pointer  to  indicate  a  course  representing  the 
sum  of  measured  track  angle  error  and  drift,  angle 
(TKE  +  DA). 

□  The  drift  angle  cursor. 

□  The  lateral  deviation  bar,  which  displays  a  linear  dis¬ 
placement  representing  cro  .s- track  error  (XTK). 

u  Ground  speed  and  distance. 

□  Engine  pressure  ratio  command. 

The  baseline  system  interconnect  block  diagram  is  shown  in  Figure  39. 

For  all  system  configurations,  the  mission  management  computer  fuel 
savings  computer  and  in-flight  data  recorder  for  flight  test  analysis  are 
mounted  in  the  equipment  bays. 
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COPILOT  INSTRUMENTS 


4.3. 2.2  Phase  Two  Configuration:  Interim  System 


For  Phase  Two,  the  system  will  be  expanded  by  addition  of  a 
multifunction  display  that  will  have  graphics  capability  for  HSI  and 
HSD/VSD  information.  The  HSD/VSD  will  show  the  horizontal/vertical 
projections  of  the  flight  plan  with  defining  waypoints  as  well  as  present 
position  of  the  aircraft  and  time  prediction. 

This  MFD  will  also  have  capability  for  overlay  of  weather  radar  signals 
on  the  HSD  information  to  show  relation  of  the  projected  flight  path  to 
local  weather  conditions. 

The  interim  system  may  be  interfaced  to  the  Speckled  Trout  to  provide 
copilot  advisories  only  for  manual  control  or,  as  an  option,  can  be  dir¬ 
ectly  interfaced  to  the  AFCS/autothrottle  to  provide  completely  automatic 
guidance  and  control  of  the  vehicle. 

The  interface  diagram  of  the  interim  system  is  shown  in  Figure  40. 

4. 3. 2. 3  Phase  Three  Configuration:  All-up  MMS 

The  final  all-up  mission  management  configuration  is  a  dual  re¬ 
dundant  system.  To  upgrade  the  interim  system  to  the  final  configura¬ 
tion  requires  the  following  additional  units. 

□  Mission  Management  Computer  (optional) 

□  Control/Display  Unit 

□  Multifunction  Display  Unit 

□  Fuel  Management  Panel 

The  additional  mission  management  computer  is  optional  since  it  is  not  a 
flight  safety  critical  item  and  may  be  required  only  as  a  backup  bus 
controller  if  the  system  elements  are  integrated  by  a  M1L-STD-1553  data 
bus.  As  an  alternate  solution,  the  FSAS  computer  function  may  be  ex¬ 
panded  to  act  as  a  backup  bus  controller. 

The  additional  CDU  and  MFD  are  to  provide  the  mission  management 
system  interface  to  the  pilot's  side  of  the  cockpit. 

The  fuel  panel  is  miniaturized  with  respect  to  the  existing  equipment. 
This  allows  space  on  the  cockpit  control  panel  for  the  additional  elec¬ 
tronic  displays.  This  fuel  panel,  however,  provides  an  additional  im¬ 
portant  function  of  controlling  fuel  out  of  the  tanks  so  that  the  c.g.  of 
the  tanker  is  maintained  at  approximately  ideal  location  for  minimum 
drag  and  hence  minimum  fuel  consumption. 
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Figure  41  is  the  interface  diagram  for  the  complete  flight  test  mission 
management  system, 

4.3.3  Flight  Test  System  Installation  Plan 

Although  the  final  installation  configuration  must  be  determined 
through  a  formal  Space  Allocation  Requirement  Program  (SARP)  and  a 
formal  Cockpit  Configuration  Review  Board  (CCRB),  an  on-site  inspec¬ 
tion  of  the  Speckled  Trout  aircraft  was  conducted  in  order  to  recom¬ 
mend  a  preliminary  plan  for  the  hardware  installation. 

4.3.3. 1  Avionics  Equipment  Bay  Installation 

The  avionics  equipment  bays  directly  aft  of  the  flight  deck,  as 
confirmed  by  Speckled  Trout  project  personnel,  is  the  appropriate  area 
for  installation  of  the  "B"  kit  component  of  the  Tanker  Mission  Manage¬ 
ment  System.  However,  because  these  avionics  equipment  bays  are 
space  critical,  and  the  space  available  for  equipment  installation  fre¬ 
quently  changes,  we  currently  have  no  specific  equipment  mounting 
recommendations.  These  specifics  will  be  established  during  the  instal¬ 
lation  design  phase  of  the  proposed  flight  test  program, 

4. 3. 3. 2  Cockpit-Mounted  Controller/Display  Units 

The  Speckled  Trout  aircraft  cockpit  was  evaluated  for  installa¬ 
tion  of  the  critical  control  and  display  equipment  required  for  each 
avionics  system  required  for  the  proposed  three-phase  flight  demonstra¬ 
tion  program,  The  modularity  of  the  flight  test  program  and  equipment 
allows  an  installation  plan  that  is  compatible  with  limited  downtime  avail¬ 
ability  of  the  aircraft  and  the  necessity  of  maintaining  the  aircraft  in 
almost  continuous  flight-ready  status. 

Phase  One  cockpit  equipment  installation  tasks  are  minimal,  requiring 
only  the  following: 

a.  Relocating  the  Collins  R-NAV  system  CDU  from  the  right  side 
to  the  left  side  of  the  center  console. 

b.  Installation  of  the  Tanker  Mission  Management  System  CDU  in 
the  space  vacated  by  relocating  the  Coiims  CDU . 

c.  Replacement  of  the  existing  EPR  indicators  with  EPRs 
equipped  with  driven  bugs. 

d.  "A"  kit  wiring  and  switching  for  interfacing  the  TMMS  with 
the  copilot's  existing  MACH/IAS  indicator  and  HS1. 

This  configuration  is  illustrated  in  Figure  43. 
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INTERFACE  DIAGRAM, 


Phase  Two  cockpit  equipment  installation  tasks  are  the  following: 

a.  Relocation  of  the  Collins  R-NAV  map  control  to  the  left  side 
of  the  center  console. 

b.  Replacement  of  the  pilot/copilot  color  weather  radar  indicator 
with  a  multifunction  display  (MFD)  and  MFD  controller. 

c.  "A"  kit  additions  to  interface  the  Digital  Flight  Guidance  Sys¬ 
tem  with  the  TMMS. 

This  configuration  is  illustrated  in  Figure  43. 

Phase  Three  involves  the  replacement  of  the  Collins  R-NAV  system  with 
a  full  dual  channel  mission  management  system,  requiring  the  following 
installation  actions  (see  Figure  44). 

a.  Replacement  of  the  copilot's  HSI  with  the  MFG  installed  for 
Phase  Two. 

b.  Replacement  of  the  pilot's  HSI  with  a  second  MFD. 

c.  Replacement  of  the  pilot's  Collins  R-NAV  CDU  with  a  second 
TMMS  CDU,  and  the  removal  of  the  R-NAV  map  display. 

Note 

For  the  Speckled  Trout  application,  it  may  be  more 
appropriate  to  install  the  second  TMMS  CDU  at  the 
navigator's  station,  along  with  a  third  MFD  and 
MFD  controller. 

d.  Removal  of  the  Collins  R-NAV  CDU,  and  RCA  Primus  color 

radar  indicator  and  controller  from  the  navigator's  station. 

e.  "A"  kit  wiring  to  interface  the  pilot's  MFD  and  CDU  with 

other  units  of  the  TMMS.  Also,  "A"  kit  switching  for 

selection  of  active  TMMS  channel. 

f.  Addition  of  LSI's  miniaturized  fuel  management  panel. 

4.4  FLIGHT  TEST  OPERATIONS 

The  modular  characteristics  of  system  mechanization,  as  defined  in 
Section  4.3,  allows  a  sequence  of  flight  tests,  as  shown  in  Figure  45, 
to  demonstrate  the  functional  and  operational  capabilities. 
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FLIGHT  TEST  SEQUENCE 
FIGURE  45 


Test  sequence  1,  using  only  the  baseline  system,  will  demonstrate  the 
basic  space/time  command  functions  of  the  system  with  the  control  loops 
for  guidance,  time  control,  and  performance  closed  through  the  pilot. 

Test  sequence  2,  with  the  upgraded  interim  system,  will  repeat  the 
basic  time/space/performance  control  flight  test  program  but  with  auto¬ 
matic  control  by  coupling  the  system  outputs  directly  to  the  autopilot/ 
autothrottle  system.  In  addition,  the  interim  system,  by  incorporation 
of  a  multifunction  display,  will  allow  flight  test  demonstration  of  in¬ 
flight  mission  redirect  capabilities. 

Test  sequence  3,  with  the  full  mission  management  system  installed,  will 
demonstrate  all  of  the  functional  and  operational  capabilities,  including 
time/fuel  trades  by  the  pilot  manager,  and  operation  without  a  navigator 
as  required  for  the  KC-135  mission. 

All  flight  tests  must  be  conducted  on  a  range  that  has  a  TACAN  station 
at  one  or  more  designated  waypoints  and  that  is  instrumented  to  track 
position  and  time  of  the  aircraft  along  selected  portions  of  the  flight 
trajectory.  The  position  and  timing  errors  of  the  ground  tracking 
system  must  be  small  compared  to  the  accuracy  requirements  of  the  avi¬ 
onics.  An  initial  estimate  of  allowable  errors  of  the  tracking  system  is 
±  500  feet  position  and  ±  1  second  time, 

4.4.1  Test  Sequence  1 

To  test  the  integrated  performance/guidance  capability  of  the 
baseline  system  to  control  aircraft  position/time,  the  flight  profile  must 
be  configured  to  exercise  the  detailed  functions.  Then  the  profiles 
must  incorporate  the  following  characteristics: 

□  VNAV 

a.  Takeoff  at  EPR  limit. 

b.  Climb  at  maximum  rate  for  noise  abatement  profile  to  in¬ 
termediate  flight  level. 

e.  Climb  at  maximum  economy  from  intermediate  to  cruise 
flight  level. 

d.  Cruise  at  maximum  economy. 

e.  Change  to  new  cruise  flight  level. 

f.  Cruise  to  top  of  descent  and  descend  on  maximum 
economy . 

g.  Appraoch  on  ILS  beam, 
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□  LNAV 


a.  Capture  initial  flight  plan  (i.e.,  initial  flight  leg  will  be 
displaced  in  direction  by  up  to  45  degrees  from  the  run¬ 
way  direction). 

b.  At  least  one  each  of  the  following  types  of  waypoints: 

□  Fixed  with  direct  flyover.  Approach  and  departure 
headings  will  differ  by  more  than  45  degrees, 

□  Fixed  without  flyover.  Approach  and  departure 
headings  will  differ  by  more  than  60  degrees. 

□  Floating  waypoint  defined  by  the  specified  flight 
path  heading  intersecting  a  VOR  radial. 

c.  Racetrack  pattern  for  rendezvous  with  receiver  aircraft. 

d.  Manual  override  of  automatic  guidance/control  with  subse¬ 
quent  return  and  recapture  of  the  flight  plan. 

e.  Typical  holding  pattern  at  a  terminal. 

f.  Approach  and  capture  of  IGS  landing  beam. 

4.4.2  Test  Sequence  2 

The  flight  test  profiles  of  sequence  2  will  be  the  same  as  for  se¬ 
quence  1,  except  that  on  at  least  some  of  the  tests,  the  pilot  will 
modify  the  flight  plan  on  the  following  schedule: 

a.  On  cue  from  data  stored  in  the  computer  as  part  of  the  flight 
plan,  the  pilot  will  investigate  alternate  LNAV  waypoints  and 
VNAV  flight  levels  to  circumvent  a  specified  threat  area  that 
is  identified  on  the  initial  flight  plan.  This  threat  could,  for 
example,  be  a  severe  weather  area.  The  criterion  for  selec¬ 
tion  of  alternate  routes  will  be  the  time  required  to  reach  the 
designated  ARCP. 

b.  On  cue  from  data  stored  on  the  computer  as  part  of  the  flight 
plan,  the  pilot  will  investigate  alternate  mission  termination 
points  to  accommodate  unexpected  shutdown  of  the  landing 
field  specified  in  the  initial  flight  plan.  The  criterion  for  al¬ 
ternate  selection  will  be  minimum  fuel  requirements  (maximum 
fuel  available  for  offloading)  to  complete  the  mission  after 
achieving  ARCP/ARCT. 
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c.  After  excursion  from  the  initial  flight  plan,  by  manual  over¬ 
ride  of  the  aircraft  control,  return  by  Direct  To  some  future 
waypoint  in  the  original  plan.  The  criterion  'or  selection  of 
the  waypoint  is  achieving  ARCP  within  ARCT  allowable  time 
deviations  and  with  minimum  fuel  expenditure. 

4.4.3  Test  Sequence  3 

The  primary  purpose  of  test  sequence  3  is  to  demonstrate  the 
total  functional  capability  of  the  complete  mission  management  system. 
This  is  in  contrast  to  sequences  1  and  2  that  are  designed  primarily  to 
test  the  functional  capability  of  the  IFTC  algorithms. 

A  secondary  purpose  of  test  sequence  3  is  to  investigate  typical  mission 
success  probability  when  operating  without  a  navigator.  The  redundant 
controller/displays  of  the  complete  mission  management  system  makes 
possible  a  flight  test  that  specifies  a  realistic  pilot/copilot  work  assign¬ 
ment,  allowing  both  crew  member's  to  interactively  operate  the  system. 

The  actual  flight  test  plan  should  simulate,  as  nearly  as  possible  but 
with  shortened  time  schedules,  the  typical  mission  scenarios  used  to 
generate  the  system  requirements  (Section  4.2). 

4.5  AIRCRAFT  INSTALLATION  AND  CHECKOUT 

Design  based  on  successful  experience  and  complete  documentation 
of  both  the  modifications  to  the  aircraft  and  the  adapter  A-Kits  assures 
straightforward  installation  of  the  avionics  system  in  the  aircraft. 

The  modification  design  will  assure  the  physical,  functional,  and  opera¬ 
tional  integration  of  the  aircraft  with  the  total  avionics  system  com¬ 
prised  of  elements  currently  on  the  aircraft  and  provided  by  the  Air 
Force,  and  those  which  LSI  designs  and  develops. 

The  total  modification  will  be  documented  in  a  Class  II  Modification 
Document,  Part  I  and  Part  II. 

4.5.1  Class  II  Modification  Document  Details 

Part  I  of  the  Modification  Document  will  contain  only  preliminary 
design  data.  Part  II  will  contain  this  design  data,  expanded  and  de¬ 
tailed,  as  well  as  pre-airworthiness  flight  test  data  and  post-airworthi¬ 
ness  flight  test  results . 
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The  design  data  includes  drawings,  analyses  and  airworthiness  flight 
test  plans.  Detailed  information  shown  on  the  drawings  includes  the 
following: 

a.  Mounting  methods,  mounting  hardware,  supports,  fixtures, 
and  mounting  position  of  each  unit. 

b.  The  specific  location  of  each  unit  in  the  aircraft. 

c.  Sway  space  clearance  provided  for  each  unit  installed. 

d.  Avionics  equipment,  airframe  material,  or  hardware  to  be  re¬ 
located  or  modified  due  to  the  installation. 

e.  List  of  materials  or  parts  list  required  for  installation,  The 
list  will  include  the  equipment  to  be  installed  and  the  installa¬ 
tion  material  or  hardware,  including  interconnecting  cables. 

The  engineering  analyses  in  the  design  data  include  hazard  or  flight 
safety,  mass  properties  (weight  and  balance),  structures,  electrical 
load,  electromagnetic  compatibility  and  stability  control. 

The  flight  test  plan  for  verification  of  airworthiness  of  the  modified  test 
vehicle  contains  test  objectives,  a  resolution  of  Air  Force  and  LSI  re¬ 
sponsibilities  in  accomplishing  the  flight  test,  test  criteria  and  proce¬ 
dures,  including  profiles  to  demonstrate  the  operational  airworthiness  of 
the  vehicle,  and  test  schedules.  Test  data  requirements,  data  reduc¬ 
tion  and  analysis  requirements,  and  flight  test  format  are  part  of  this 
test  plan. 

This  test  plan  will  be  preliminary  in  Part  1  and  updated  and  modified 
for  the  final  version  in  Part  II  of  the  Class  II  Modification  Document. 

In  addition  to  updated  versions  of  major  sections  of  Part  1,  the  Part  II 
Modification  Document  also  contains  Operating/Maintenance  Instructions, 
lists  of  applicable  supporting  documentation  and  spare  parts,  pre-air- 
worthiness  flight  test  data  and  post-airworthiness  flight  test  results. 

A  partial  flight  manual  will  provide  the  necessary  ground  and  flight 
operation,  maintenance,  and  inspection  instructions  for  the  modified  test 
aircraft . 

4.5.2  Aircraft  Installation  Design 

The  aircraft  Class  II  modification  to  install  and  integrate  the  avi¬ 
onics  system  into  the  test  aircraft  must  be  in  accordance  with 
MIL-P-27733  and  safety  of  flight  requirements. 
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The  pilot's  and  copilot's  main  instrument  panels,  as  shown  in  Figure  44, 
will  be  redesigned.  All  new  Controller/Display  equipment  will  be  sup¬ 
plied  as  GFE.  New  instrument  panels  will  be  provided  as  part  of  the 
A-Kit.  Electrical  interface  of  this  equipment,  as  well  as  all  other  ele¬ 
ments  of  the  total  avionics  system ,  will  be  included  in  our  integration 
design . 

Flight  safety  during  the  program  is  assured  by  the  modularity  of  flight 
test  avionics  and  test  plan  and  by  an  engineer  specifically  assigned  the 
responsibility  for  coordinating  all  safety-related  activities,  analyses, 
and  reports.  The  safety  efforts  must  begin  at  the  start  of  the  program 
and  carry  through  completion.  The  efforts,  results,  and  conclusions 
will  be  documented  in  a  Preliminary  Hazard  Analysis  and  a  Final  System 
Hazard  Analysis. 

4.6  FLIGHT  TEST  SUPPORT 

LSI  will  provide  support  to  the  Air  Force  for  modification  of  the 
test  aircraft  and  installation  of  all  Group  A  and  Group  B  equipment  as 
defined  by  the  Class  II  Modification  Document.  Specific  support  tasks 
(Figure  46)  include  ground  testing,  systems  calibration  and  checkout, 
preparation  of  the  aircraft  for  functional  check  flights,  actual  flight 
testing,  integration  support  of  advanced  systems,  equipment  mainten¬ 
ance,  environmental  assessment  and  field  support. 

SPECKLED 


fflTTTT' 


SYSTEM  DEMO  TASK  FLOWCHART  -  PHASE  111 

FIGURE  46 
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Air  Force  personnel  at  Detachment  1,  4950th  Test  Wing  (Speckled 

Trout)  will  actually  accomplish  the  modification,  installing  both  the 
equipment  which  LSI  will  supply  and  equipment  furnished  by  the  Air 
Force, 

4.6.1  Aircraft  Modifications 

LSI  support  will  assure  that  aircraft  modifications  comply  with  all 
applicable  military  specifications,  the  approved  Class  II  Modification 
Document,  and  Air  Force  Quality  Control  standards  and  procedures. 

The  applicable  military  specifications  are  MIL-P-27733,  MIL-1-45208, 
MIL-STD-882,  and  T.O.  1-1A-14.  The  Class  II  Modification  Document 
will  be  modified  as  necessary  as  the  airframe  modification  progresses  to 
incorporate  the  actual  physical  configuration  of  the  system. 

The  Air  Force  will  prepare  the  test  aircraft  and  remove  any  equipment 
from  the  test  aircraft  that  is  necessary  to  accomplish  installation  of  the 
respective  Group  A  and  Group  B  kits. 

All  Group  A  components  of  the  system  will  be  installed  in  accordance 
with  AF-appro'  'd  Class  II  modification  documentation,  MIL-I-45208,  and 
Air  Force  Quality  Assurance  procedures.  LSI  and  government,  inspec¬ 
tions  will  be  conducted  during  modification. 

4.6.2  Ground  Checkout 

A  careful  ground  checkout  --  including  EMI  tests,  system  cali¬ 
bration  and  functional  tests  --  assures  that  the  test  aircraft  is  ready 
for  the  Air  Force  functional  test  flight. 

EMI  testing  verifies  that  EMI  integrity  within  the  aircraft  is  not  affected 
by  the  equipment  installation. 

System  functional  testing  verifies  proper  installation,  alignment,  calibra¬ 
tion,  and  operation  of  all  systems  and  disturbed  systems, 

These  final  ground  flight  preparation  checks  include  weight  and  balance 
check,  compass  rose  swing  (if  required),  hydraulics  check,  landing 
gear  check,  90-day  calendar  inspection  (if  required),  brake  and  tire 
check,  nose  wheel  steering  check,  and  fuel  level. 

The  test  aircraft  is  now  ready  for  a  functional  check  flight  by  Air 
Force  flight  crew  to  verify  system  operation  prior  to  formal  flight  test¬ 
ing.  LSI  will  correct  any  deficiencies  found  during  this  functional 
check  flight.  This  correction  of  deficiencies,  if  any  are  found,  and  an 
acceptable  functional  flight  test  complete  the  modification/installation 
phase . 
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4.6.3  Data  Analysis 


LSI  will  develop  test  criteria  and  provide  data  analysis  support 
that  will  maximize  the  demonstration  results  with  a  minimum  number  of 
flight  hours, 

For  this  task,  on-site  engineers  and  technicians  will  be  established, 
supported  by  engineers  in  Grand  Rapids  as  required,  to  assure  effec¬ 
tive  operation  of  the  avionics  system. 

Flight  test  criteria  will  be  developed  to  demonstrate  the  system  airworth¬ 
iness  .  This  task  includes  definition  of  all  the  operational  data  require¬ 
ments,  including  source,  function,  and  analysis  methods. 

The  basic  performance  criteria  are  consistent  non-glitching  operation  of 
the  integrated  system  and  consistent  output  of  the  data  bus  parameters 
to  the  Data  Recording  System. 

These  criteria  are  incorporated  as  the  flight  test  plan  of  the  Class  II 
Modification  Document. 

Following  the  flights,  LSI  will  evaluate  the  data  obtained  by  the  on¬ 
board  recording  equipment  as  well  as  the  observations  of  the  flight 
crew.  This  flight  test  data  will  be  summarized  for  inclusion  in  both  the 
interim  and  final  flight  test  reports. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  shows  that  the  complexity  of  generating  a  four-dimensional 
trajectory  using  optimal  techniques  was  too  great  to  allow  more  than  a 
cursory  examination  of  the  problem  within  the  time  and  money  con¬ 
straints  of  the  program.  The  principal  problem  was  in  defining  an  opti¬ 
mization  algorithm  that  could  readily  satisfy  the  initially  defined  cost 
function . 

A  two-step  approach  of  separating  the  horizontal  plane  and  vertical 
plane  projections  of  the  desired  trajectory  was  investigated  briefly . 
This  approach,  while  analytically  more  amenable,  was  not  attempted  be¬ 
cause  of  the  inherent  reductions  in  performance. 

Optimal  control  of  the  aircraft  by  a  linear  quadratic  regulator  to  the  tra¬ 
jectory  generated  by  classical  algorithms  was  investigated.  Although 
technically  feasible,  this  technique  proved  to  be  too  costly  in  terms  of 
on-board  avionics  computer  capabilities  to  justify  the  probable  degree  of 
performance  improvement. 

A  three-phase  plan  for  flight  verification  of  the  classical  Flight  Trajec¬ 
tory  Control  algorithms  and  for  demonstration  of  their  operational  capa¬ 
bilities  was  defined.  This  flight  test  plan  will  utilize  the  Speckled 
Trout  test  aircraft  and  existing  avionics  hardware  with  modified  soft¬ 
ware  . 

This  plan  is  recommended  as  a  program  follow-on  because  it  will  do  the 
following: 

a.  Provide  an  experimental  data  base  for  the  utility  of  the 
Flight  Trajectory  Control  concept  --  particularly  in  moni¬ 
toring  preplanned  time  on  target  with  in-flight  redirect; 

b.  Establish  the  computer  resources  of  the  on-board  avi¬ 
onics  to  implement  the  Flight  Trajectory  Control  concept; 

c.  Establish  operational  characteristics  of  the  concept  for 
application  to  an  integration  with  other  developmental  ca¬ 
pabilities  of  advanced  aircraft  such  as  terrain-following 
and  integrated  fire/flight  control. 
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APPENDIX  A 


AIRCRAFT  TRANSFER  FUNCTIONS 


For  both  flight  conditions,  the  values  of  the  aerodynamic  coefficients 
are  the  following: 


Aerodynamic  Coefficient 

Flight  Condition 

Low  Q 

High  Q 

C 

m6 

-.688  (1/rad) 

-.481  (1/rad) 

e 

CT 

L6 

.223  (i/rad) 

.149  (1/rad) 

e 

C 

m* 

-4.7  (1/rad) 

-5.4  (1/rad) 

a 

C 

ma 

-1.09  (1/rad) 

-1.32  (1/rad) 

CL 

4.56  (1/rad) 

5.54  (1/rad) 

a 

CD 

.049 

.0152 

C 

mq 

-15.1  (1/rad) 

-15.25  (1/rad) 

Solving  for  the  stability  derivative,  the  following  results  are  achieved: 


Stability  Derivative 

Flight  Condition 

Low  Q 

High  Q 

-.93 

-2.61 

e 

‘  z6 

-10.23 

-27.4 

e 

-4 

-8.48x10 

-5.05x10 

-.00536 

-.00938 

zw 

-.769 

-.927 

M 

-.749 

-1.092 

q 

Therefore,  for  the  low  Q  flight  condition, 

,  v  _  -.921(8+. 717) 

6  (S)  ~  +  fT75s  +  2.05 


and  for  the  high  Q  flight  condition, 

,  v  _  2,6(s+.83) _ 

6'  ~  bz  +  2.4s  +  8.17 

e 

10.2  ~  (s)  TRANSFER  FUNCTION 


The  3-degree-of-freedom  lateral- directional  equations  of  motion  are 
the  following: 


The  lateral- directional  aerodynamic  coefficients  for  both  flight  conditions 
are  given  below. 


Aerodynamic  Coefficient 

Flight  Condition 

Low  Q 

High  Q 

cn 

-.026  (1/rad) 

.021  (1/rad) 

P 

C 

ye 

-.676  (1/rad) 

-.756  (1/rad) 

S 

-.186  (1/rad) 

-.17  (1/rad) 

p 

c* 

-.364  (1/rad) 

-.306  (1/rad) 

p 

C£ 

.203  (1/rad) 

.129  (1/rad) 

r 

C 

QP 

.128  (1/rad) 

.142  (1/rad) 

C 

-.158  (1/rad) 

-.148  (1/rad) 

r 

C 

.211  (1/rad) 

.179  (1/rad) 

% 

6r 

.0099  (1/rad) 

.0184  (1/rad) 

% 

a 

.0253  (1/rad) 

.0198  (1/rad) 

C 

n, 

6r 

-.0843  (1/rad) 

-.0733  (1/rad) 

c. 

nA 

6a 

0  (1/rad) 

.0034  (1/rad) 

Solving  for  the  stability  derivatives,  the  following  results  are  achieved: 


Stability  Derivative 

Flight  Condition 

Low  Q 

High  Q 

N 

P 

-.0318 

.0386 

Y 

V 

-.0269 

-.0156 

-2.31 

-8.46 

L 

P 

-1.07 

-1.305 

L 

r 

.6 

.55 

NP 

.684 

3.04 

N 

r 

-.201 

-.272 

% 

.0352 

.0431 

.  123 

.916 

•  6a 

.314 

.985 

N, 

6r 

-.45 

-1.57 

N, 

6a 

0 

.073 

| 
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Substituting  the  stability  derivatives  into  the  lateral- directional  equations 
of  motion,  we  have  the  following  for  the  low  Q  flight  condition: 


S+.0269 

2.31 

-.684 


-.117 
•2+l .07s 
- .0319*+. 0318s 


and  for  the  high  Q  flight  condition: 


s+,0156 

8.46 

-3.04 


The  aircraft  transfer  function  t-(s)  is  obtained  using  the  IBM  370  program 
LISA.  °R 


For  the  low  Q  flight  condition: 

t  /  \  _  -.447(s+l. 25) [(s-.0934)z+(.437)2] 

6r(s)  "  (s+1 . 26)(s+. 00589) [ (s+ . 00856) 2+( . 923 


and  for  the  high  Q  flight  condition: 


-1 .54(s+l .44) [ (s- . 0976) 2+ ( . 436) 2 ] 
Js+1.39Ks+.  00668)  [(8+.0913)2+(1. 69)2] 


10.3  A  (s)  TRANSFER  FUNCTION 

a  r+6R 

The  equation  relating  r(s)  to  ^(s)  is  the  following: 

Vs)  "  (  Q  ^^s+lOs+l  ^S+15 . 7  r(s) 

Modifying  the  lateral- directional  transfer  function  matrix,  we  have  the 
following  for  the  low  Q  flight  condition: 
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s+,0269  -.117 

2.31  s2+1.07s 


1.  -.0352 

-.074s-. 6  -.123 


-.684  -.0319s2+. 0318s  S+.201  .45  rCs) 

0  0  -(31 ,400/Q)s  s3+26.7s2+182.7s+157j[_6R(s) 

and  for  the  high  Q  flight  condition: 


6  (s) 
A 


-.0422 

1 . 

-.0431 

s2+l . 035s 

-.074s-. 55 

-.916 

-.0319s2-. 0386s 

S+.272 

1.57 

0 

(-31  ,400/Q)s 

s3+26.7s2+182 

Vs> 


The  aircraft  transfer  function,  r-^(s)  >  is  obtained  by  using  the 

°A  ^ 

analysis  program  LISA. 


a  r+6 


For  the  low  Q  flight  condition: 

_ .  315  ljs+  .  l46)2+( •  547 ) 2  ]  (s+1 . 75)  (s+6  ■ 87) (s+1 7 . _1 ) _  . 

r*6  ~  [ (s+ .0468)2+( . 7 13) 2 ] (s+1 . 25) (s+ .00463) (s+6 .  29) (s+3 . 09) (s+1 7 . 25) 


Am 


and  for  the  high  Q  flight  condition: 

•  9 ? 3 1  (s+  ■  628) 2 +  (1  .54) 2  1  (s+1  .  75)  (s+6 . 87)  (s+1 7  .  1) _ 

r->6  ~  ( (s+. 395)2+( 1.35)2] (s+1. 61 )(s+. 00556) (s+1 . 75) (s+7 . 03) (s+17 . 1) 


APPENDIX  B 

STOCHASTIC  THREAT  MODEL 


I 


•i 


fi 
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STOCHASTIC  THREAT  MODEL 


In  order  to  develop  a  robust  algorithm  for  resource  allocation,  it  is 
necessary  to  take  into  account  the  stochastic  nature  of  the  threat  en¬ 
vironment.  If  the  solution  were  to  consider  only  those  threats  which 
absolutely  were  known  to  exist,  then  it  would  leave  the  aircraft  vulner¬ 
able  to  "pop-up"  threats.  On  the  other  hand,  if  the  algorithm  were  to 
consider  all  possible  pop-up  threats  which  conceivably  might  appear, 
then  it  would  tend  to  avoid  threats  which  in  fact  never  materialize. 
Such  a  worst-case  policy  is  also  inappropriate.  What  is  required  is  an 
analytical  procedure  which  results  in  a  proper  trade-off  of  flight  path 
between  known  present  threats  and  possible  future  threats.  Such  a 
procedure  is  developed  in  this  section. 

Three  types  of  information  are  assumed  to  be  available  to  the  the  air¬ 
craft  for  each  known  present  and  potential  future  threat.  They  are 
as  follows: 

1.  a  probability  that  the  threat  actually  exists; 

2.  a  probability  distribution  for  the  location  of  the  threat 
given  that  it  does  exist;  and 

3.  a  threat  model  specified  by  the  conditional  probability  of 
aircraft  survival  given  that  the  threat  actually  exists 
and  is  at  a  known  location,  as  a  function  of  the  air¬ 
craft's  position  and  velocity  relative  to  the  threat. 

The  equations  are  derived  which  transform  this  information  into  an  un¬ 
conditional  threat  model.  Then,  certain  simplifying  assumptions  are 
described  which  allow  the  rapid  real-time  computation  of  the  uncondi¬ 
tional  threat  model. 

20.1  DERIVATION  OF  GENERAL  EQUATIONS 

For  each  threat,  assume  that  a  value,  for  the  probability  that  the 
threat  exists,  P(X),  is  known.  P(X)  will  be  near  1  for  presently 
known  threats,  and  will  be  less  for  postulated,  future  pop-up-type 
threats.  Furthermore,  a  probability  distribution  of  each  threat's  loca¬ 
tion  is  assumed  to  be  known.  The  form  of  the  distribution  is  arbi¬ 
trary,  although  simple  analytic  expressions  such  as  uniform  or  Gaussian 
distributions  are  reasonable  assumptions  and  certainly  more  tractable 
than  quantized  representations  requiring  numerical  integration.  This 
distribution  is  expressed  conveniently  in  polar  coordinates  (see  Figure 
B-l),  relative  to  the  man  threat  location  denoted  f(R,0).  The  threat 
model  is  a  function  of  three  position/velocity  variables,  as  follows: 
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AIRCRAFT 

VELOCITY 


P(S|p,M,A,X) 


Probability  of  aircraft  survival  per  unit 
time  for  horizontal  separation  p,  aircraft 
approach  angle  p,  aircraft  altitude  z, 
resource  allocation  A,  and  conditioned 
upon  the  threat  existing. 

The  desired  threat  model  which  is  to  be  derived  below  is  the  uncondi¬ 
tional  model  (i.e.,  no  dependence  on  the  random  variable  X)  with  re¬ 
spect  to  t'\e  mean  of  the  threat  distribution: 

P(S|r,y,z,A)  =  Probability  of  aircraft  survival  per  unit 

time  for  horizontal  separation  between  the 
aircraft  and  the  mean  of  the  threat  dis¬ 
tribution  r,  angle  of  aircraft  velocity  with 
respect  to  the  line  joining  the  aircrait  and 
the  threat  center  y,  aircraft  altitude  z 
(see  Figure  1),  and  resource  allocation  A. 

The  first  step  in  the  derivation  is  to  distribute  the  unconditional  prob¬ 
ability  over  the  two  mutually  exclusive  events  of  the  threat's  existence 
and  its  nonexistence  P(X). 

P(S|r,y,z,A)  =  P(S|r,y,z,A,X)P(X)  +  P(S| r , y,z,A,X)P(X) . 

But,  if  the  threat  does  not  exist,  then  the  survival  probability  must  be 
unity.  So: 


P(S|r,y,z,A) 

P(S|r,Y>z,A,X)P(X)  +  1  -  PCX). 

(B-l) 

The  next  step  is  to  distribute  the  conditional  probability  appearing  on 
the  right  side  of  the  above  equation  over  all  possible  locations  for  the 
threat,  as: 

P(S | r , y,z ,A,X) 

- 

JJf(R,e)P[S|p(R,0),p(R,0),Z,A,X]dRd0 

(B-2) 

where  (using  the  laws  of 
p(R,0) 

sines 

and  cosines  and  Figure  B-l) 

R2  +  rz  -  2Rrcos0 

(B-3a ) 

P(R,e) 

= 

Y  +  «KR,0) 

(B-3b) 

where 

sinijj 

= 

Rsin0 

P 

(B-3c) 

cos 

zz 

r-Rcos0 

n 

(B3-d) 
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Thus,  a  general  solution  to  the  problem  exists.  Given  the  three  items 
of  information  about  a  threat,  the  following  procedure  may  be  used  to 
develop  the  unconditional  threat  model.  For  a  given  value  of  r,Y,z 
and  A,  the  conditional  probability  P(S|r,Y,z,A,X)  may  be  calcu¬ 
lated  by  numerical  integration  of  equation  B-2,  using  equations  B-3a 
through  B-3d.  Then  the  desired  unconditional  probability  may  be  ob¬ 
tained  by  using  equation  B-l. 

20.2  SIMPLIFICATION  OF  GENERAL  EQUATION 

In  the  preceding  section,  relationships  between  the  conditional  and 
unconditional  probabilities  of  survival  were  derived  in  general.  In  this 
section,  the  general  equations  will  be  simplified  so  that  rapid  real-time 
computation  is  possible. 

The  first  thing  to  note  is  that  the  mission-planning  dynamic¬ 
programming  algorithm  does  not  operate  on  the  probability  of  survival 
directly,  but  rather  on  the  negative  logarithm  of  the  survival  probabil¬ 
ity.  Therefore,  the  following  notational  convention  is  adopted: 

E(variables)  -  -£n  (P(S | variables) ]  (B-4) 

The  quantized  unit  of  time  proposed  for  the  mission  planner  will  be  a 
short  enough  length  of  time  that  all  survival  probabilities  are  nearly  1. 

Hence,  the  following  approximation  is  of  excellent  numerical  accuracy. 

E(variables)  =  1  -  P(S | variables)  (B-5) 


Applying  the  above  relationship  to  equation  B-l, 


P(S I r , y,A,z) 

1  -  P(S | r ,y,A,z) 
E(r,Y,A,z) 

and  then  to  equation  B-2, 
E(R,y,z,A,X) 


E(r,Y,z,A,X) 


1  -  P(X)  [ 1 -P (S | r , Y>A,z,X) 1 
P(X)  [1  -  P(S|r,Y,z,A,X)] 

P(X)  E(R,y,z,A,X) 

1  -  j/f(R,e)p(s|p,p,z,A,x)dRde 
JJf(R,e)dRde  -  J/f (R,8)P(S|p,0,z,A,X)dRd0 
JJf(R,e)[l  -  P(S|p,0,z,A,X)] dP.de 
J/f(R,e)  E(p,p,z,A,X)dRd0 


(B-6) 


(B-7) 
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and  combining  the  results  (equations  B-6  and  B-7)  give  the  following: 

E(r,Y,z,A)  S  P(X)  .f/f(R,0)E(p,p,2,A,X)dRde 

Equation  B-8  may  be  used  for  numerical  integration,  along  with  equa¬ 
tions  B-3a  through  B-3d  and  using  a  uniform  or  Gaussian  probability 
density  function  for  f(R,6)  as  appropriate.  This  procedure  is  perfectly 
general  as  long  as  the  unit  of  time  is  taken  to  be  small  enough  so  that 
the  approximation  in  equation  B-5  is  valid. 
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